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Abstract. We compute the static potential in the Gribov-Zwanziger Lagrangian as a function 
of the Gribov mass, 7, in the MS scheme in the Landau gauge at one loop. The usual gauge 
independent one loop perturbative static potential is recovered in the limit as 7 — )• 0. By contrast 
the Gribov-Zwanziger static potential contains the term 7^/(p^)^- However, the linearly rising 
potential in coordinate space as a function of the radial variable r does not emerge due to a 
compensating behaviour as r — )• 00. Though in the short distance limit a dipole behaviour is 
present. We also demonstrate enhancement in the propagator of the bosonic localizing Zwanziger 
ghost field when the one loop Gribov gap equation is satisfied. The explicit form of the one loop 
gap equation for the Gribov mass parameter is also computed in the MOM scheme and the zero 
momentum value of the renormalization group invariant effective coupling constant is shown to 
be the same value as that in the MS scheme. 
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1 Introduction. 



One of the main outstanding problems in quantum field theory is to understand the confinement 
of the quarks and gluons associated with the strong nuclear force. The underlying field theory is 
Yang-Mills or Quantum Chromodynamics (QCD) which is a non-abelian gauge theory based on 
the SU{3) colour group. Whilst its ultraviolet properties such as asymptotic freedom, [lii2j, are 
relatively straightforward to extract from perturbative calculations, it is the infrared properties 
which are difficult to truly access. These lie at the heart of the confinement problem. In early 
analyses of trying to understand how quarks and gluons are confined, it was recognised that the 
force separating these fundamental quanta of QCD behaved as a constant at large distances. In 
other words the potential energy at a colour charge separation distance of r was proportional to 
r. Known as the linearly rising potential it has been studied using lattice regularization which 
goes beyond perturbation theory. Indeed there is widespread acceptance that such numerical 
analyses do produce a linearly rising potential. See, for example, [Hill [5]. Another aspect of 
early attempts to understand such a potential has been clearly given in Mandelstam's review, 
[6], which emphasises the infrared behaviour of the gluon field, as well as in the work of others, 
[71[8l|9]. Briefly, in the ultraviolet the gluon behaves as a standard fundamental massless particle 
with a simple propagator proportional to l/p^ where p is the gluon momentum. However, such 
a behaviour can really only be valid asymptotically at high energy. Unlike the photon and 
overlooking for the moment the fact that the gluon is not a gauge invariant concept, if the gluon 
propagator's behaviour was of this form for all momenta then it could naively but erroneously 
be identified as a real physical object contrary to its accepted confined property. Instead due to 
the non-linear nature of the non-abelian gauge theory the gluon's propagator must be modified 
non-perturbatively or dynamically to a non-fundamental form without a simple pole. Given 
this, one natural form which it might take in order to fit with the linear potential is a dipole. 
The justification for this is that considering the force between coloured sources, the 
exchange of such a gluon at low energies will naturally produce a linear potential in coordinate 
space after taking a simple Fourier transform. However, whilst such a behaviour is an ultimate 
goal, the evidence for a dipole behaviour in the gluon propagator never actually manifested 
itself in subsequent studies. Indeed with the use of more powerful computers, lattice gauge 
theories have been able to probe the infrared behaviour of the gluon propagator more deeply 
at low energy. Without getting drawn too deeply for the moment into the current debate over 
which of the scaling, [inillllliaiiailllllllllllTllIHllIllSOlEl^ or decoupling solutions, 
[22l[23ll2l[25l[26l|271[28l[l9], is correct, it is fair to say that the gluon does not appear to 
have a dipole form at low momenta. Instead the propagator freezes either to zero or a finite 
non-zero value for each of these respective current scenarios. Therefore, the extraction of the 
dipole behaviour is currently in abeyance. 

The next key breakthrough in the pursuit of gluon confinement effectively dawned with 
Gribov's seminal work noting the ambiguity in gauge fixing in the Landau gauge, pTOj. In 
essence gauge fixing locally in a non-abelian gauge theory can be achieved uniquely but globally 
one will always end up with gauge or Gribov copies hindering the process. This was resolved 
by Gribov, |10) . by restricting the path integral measure to the first Gribov region, O, defined 
as that region encompassing the trivial gauge solution, = 0, where AJJ is the gluon field or 
gauge potential, such that the Faddeev-Popov operator has strictly positive eigenvalues. This 
restriction of the path integral produces a radically different gluon propagator which is clearly 
non-fundamental and vanishes in the infrared limit. Moreover, a mass parameter is introduced, 
known as the Gribov mass and denoted by 7, which is not an independent parameter of the 
theory. Indeed 7 satisfies a gap equation which can be evaluated at one loop, [TO]. One major 
consequence is that the propagator of the Faddeev-Popov ghost, which is required for the Landau 
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gauge fixing, also has a modification. Instead of a simple l/p^ form, it behaves as as 
— )• 0. Whilst this is a dipole, a Faddeev-Popov ghost clearly cannot be exchanged between 
colour sources as it is Grassmann in nature as it is required to restore unitarity. Subsequent 
progress in this area was via a series of articles in the main by Zwanziger and collaborators, |1H 
[laillllllllSllTBlinillHllTSllSOlEl]. Briefly, the drawback of the Gribov semi-classical analysis 
is that the resultant effective Lagrangian is non-local and therefore essentially inadequate for 
practical quantum computations. However, in a series of articles Zwanziger managed to localize 
the non-locality at the expense of introducing extra localizing fields, |12i [T3l [T71 I18j . The 
resulting Lagrangian was renormalizable, |18 | 1301 [3T] . meaning that it was possible to perform 
calculations with it. Indeed its renormalization is such that the known ultraviolet structure 
of Yang-Mills or QCD is totally unaltered by the extra localizing fields. Instead they only 
effectively play a role in the infrared. Subsequently, the two loop MS extension to Gribov's gap 
equation was determined in [32j verifying the Faddeev-Popov ghost enhancement at two loops 
and gluon suppression at one loop, [33]. The former property is not inconsistent with the Kugo- 
Ojima confinement criterion, [Ml EH], derived for a BRST invariant Yang- Mills theory. Indeed 
it was shown that the Grassmann localizing ghost also was enhanced at two loops, [33] . giving 
a propagator of dipole form but it is equally not directly relevant for confinement of coloured 
sources for the same reason as the Faddeev-Popov ghost. 

Given the potential for the Gribov-Zwanziger Lagrangian to be a candidate for understanding 
confined gluons we now indicate the primary aim of this article. The explicit calculations of 
[32| [33] at one and two loops have opened the possibility of calculating the potential between 
two coloured sources using the Gribov-Zwanziger localized Lagrangian, [T2l |13l [l7l [18] . One aim 
is to see whether the behaviour of the potential is significantly different from that computed in 
the usual version of QCD. Indeed this was first discussed by Susskind in [36] and other authors 
at around the same time, [371 ESj- They considered the Wilson loop definition and examined 
the energy between two coloured sources fixed in time but at a spatial separation. This is 
known as the static limit. An advantage of the Wilson loop is that it is gauge invariant and 
in principle one can extract the resulting potential in any gauge. Originally this was achieved 
in the Feynman gauge at one loop in the MS scheme, [36l [371 EH]- Subsequently, the two loop 
static potential was determined in the same gauge in [391 [40] . This calculation was subsummed 
into a more general computation by Schroder in |41[ I42j . There the full two loop potential was 
constructed in an arbitrary linear covariant gauge in ordinary (massless) perturbation theory. It 
clearly demonstrated the explicit cancellation of the linear gauge fixing parameter (and en route 
corrected a minor error in the original expression given in |39[I40]). The two loop static potential 
for a general colour state emerged later in [53| . More recently the three loop static potential has 
been the subject of interest with both the quark, |44[ I45[ [^. and purely gluonic contributions 
now available, [471 148] . Since such an extensive formalism already exists (and is comprehensively 
reviewed in [49] ) it is a rather straightforward exercise to apply it to the Gribov-Zwanziger 
Lagrangian at one loop. At the very least any perturbative potential, which will depend on 7, is 
needed since it would have to match onto the behaviour of the potential beyond the perturbative 
approximation. At this point it is worth noting that Zwanziger considered the Wilson loop in 
the Gribov-Zwanziger context using a lattice approach in [T9]. There it was argued that a linear 
rising potential could emerge if there was ghost enhancement and gluon suppression, including 
gluon propagator freezing to a non-zero value. Moreover, the string tension was proportional to 
a combination of the zero momentum values of the gluon and Faddeev-Popov ghost form factors. 

Although the static potential forms a main part of our article, we will also consider the 
structure of all the contributing fields to see whether a dipole exchange between coloured sources 
could somehow emerge in the Gribov-Zwanziger Lagrangian. In making this previous statement 
we have been careful in our wording and not mentioned the gluon. The reason for this is simple. 
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The premise that it is solely the gluon field itself which is responsible for confinement in the 
Gribov-Zwanziger context may need to be refined. This is primarily due to the appearance of the 
Zwanziger localizing ghosts in the work of [121 [T3l [T71 [T8] . These play no role in the ultraviolet 
dynamics of QCD but do become important in the infrared. There is a bosonic Zwanziger 
ghost which is spin-1 and carries non-abelian colour charge. In some ways it could be regarded 
as a gluon component since its equation of motion implies it is a non-local projection of A^. 
However, in the full Gribov-Zwanziger Lagrangian we regard it as a separate entity. Therefore, 
in principle the appearance of a dipole behaviour in the infrared between coloured sources could 
also derive from other spin-1 fields aside from the gluon itself. Indeed this is indicated very 
strongly in Zwanziger 's recent article, |50j . where Dyson Schwinger techniques are applied to the 
bosonic localizing ghost to produce an enhancement akin to that of both Grassmann ghost fields. 
Though the enhancement is actually beyond dipole being l/(p^)^ as — )• 0. However, since it 
is non- Grassmann, [50], it is evidently a much better candidate for an exchange field between 
coloured sources. Therefore, in the current context of the full one loop MS static potential in 
the Gribov-Zwanziger Lagrangian we will speculate how it might be possible to go beyond the 
one loop potential and qualitatively produce a linearly rising potential. En route we will record 
some new properties of the underlying Gribov gap equation. As emphasised by Zwanziger, in 
the Gribov-Zwanziger context the theory can only be regarded as a gauge theory when the 
Gribov mass, 7, explicitly satisfies this equation. In the original calculations of Gribov, [10], to 
produce the Faddeev-Popov ghost enhancement the main lesson was that the gap equation was 
central to seeing towards the infrared. To a lesser extent it also played a role in qualitatively 
producing the freezing to a non-zero value, [33], of a renormalization group invariant definition 
of the coupling constant based on the Landau gauge properties of the gluon ghost vertex, [5T] . 
Therefore, we believe that the gap equation, which is derived from the condition defining the 
first Gribov region, should be a central requirement to producing the infrared property of a 
linearly rising potential. In our calculations we will follow the more recent formulation of the 
Gribov-Zwanziger Lagrangian of [50] where the localizing bosonic ghosts were resolved into their 
real and imaginary parts. Whilst this will be different from the earlier loop computations of 
[32I [33] . it transpires that those results are not correct since the propagator of the real part 
of the bosonic localizing ghost was treated erroneously. Therefore, throughout the discussion 
we will re-address some of the results of [32l [33] and provide the correct details. It turns 
out that the main features such as one loop gluon suppression, two loop Faddeev-Popov ghost 
enhancement and the freezing of a renormalization group invariant effective coupling constant 
to a finite non-zero value remain unaltered. It is only the actual numerical details which are 
revised since the omission affects the finite parts of Feynman graphs upon which there are no 
non-trivial independent checks. Finally, given the current interest in the alternative infrared 
structure known as the decoupling solution, [22l [23l 1231 [25l [26l [271 [28l [29] , we wih briefiy discuss 
its status within the static potential approach. Indeed in [52] a novel way of producing the 
decoupling solution within the Gribov-Zwanziger Lagrangian was considered and later tested in 
explicit calculations, |53) . 

The paper is organised as follows. In section 2 we review the static potential formalism of [36[ 
[37 1 138 1 l39 l I40 [ HT l H2] as well as the necessary ingredients from the Gribov-Zwanziger Lagrangian, 
[I2I [T3l [TTI [TS] . Section 3 is devoted to the formal construction of the one loop corrections to 
the propagators in the reformulation of the Gribov-Zwanziger Lagrangian in terms of the real 
and imaginary parts of the localizing bosonic ghost discussed in [JO]. The explicit one loop MS 
static potential is constructed in section 4 prior to considering several new calculational aspects 
of the Gribov gap equation for 7 in section 5 including the correct two loop MS expression. The 
role the original gap equation plays in the one loop enhancement of the bosonic localizing ghosts 
is discussed in section 6 and its implications for trying to extract the linearly rising potential 
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qualitatively in the Gribov-Zwanziger static potential are also considered. Section 7 focuses on 
considering power corrections to the potential and other quantities which can be computed in 
the Gribov-Zwanziger Lagrangian as well as elementary ideas on the underlying properties of 
higher order Feynman graphs. We give our conclusions in section 8. Several appendices are 
provided. The first details the decomposition of products of colour group structure functions 
into a standard basis which is required in the discussion of the bosonic ghost enhancement. The 
next two appendices provide the explicit exact one loop corrections respectively to the transverse 
and longitudinal parts of all the 2-point functions of the spin-1 fields in the Gribov-Zwanziger 
Lagrangian. The longitudinal parts are presented as they are relevant to the current debate 
on the BRST structure of the Gribov-Zwanziger Lagrangian. The final appendix records the 
explicit values of scalar amplitudes in the decomposition of the correlation function of a Lorentz 
tensor operator involving the field strength. It extends the power correction computation of the 
analogous scalar operator considered in section 7. 



2 Formalism. 



We begin by briefiy summarizing the static potential formalism in both the original QCD and 
Gribov-Zwanziger contexts. The key is the definition of the Wilson loop where the loop is taken 
to be a rectangle. The spatial extent is length r and the time interval is denoted by T. When 
the temporal side of the rectangle is very much in excess of the spatial extent, T ^ r, then the 
potential, V{r), between two heavy static colour sources is given by 



V{r) 



lim — In ( 

T-s>oo iT 



Tr V exp ig 



fj. 



(2.11 



Here V denotes the path ordering prescription, is the gluon field and are the colour group 
generators satisfying 

[T^r''] = if^^T^ (2.2) 

where f""^^ are the structure constants. Throughout our discussion on the definition and prop- 
erties of the static potential we make the same general assumptions as have been discussed 
extensively in [Ml [37l [Ml EHl HQI HU |42] such as gauge invariance of V{r). Also, for instance, it 
is accepted that the potential of (|2.ip can be reformulated in terms of a functional integral with 
an external colour source, J^{x). In other words 



Vir) 



1 trZ[J] 
lim — — — — r 

T^oo iT trZ[0] 



where formally 



Z[J] 



VAnVipV'4)VcVc exp 



j (fx (l 



+ 



(2.3) 



(2.4) 



and the Lagrangian L will be discussed later in detail but will either be the original QCD 
Lagrangian or that of the Gribov-Zwanziger theory. The denominator, Z[0], is included as 
a normalization. The other main assumption, which also occurs in the usual QCD case and 
has been discussed in [Ml [37l [Ml [Ml HOI E], is that of the exponentiation of the values of 
the actual Feynman diagrams in order to properly extract the static potential implicit in the 
definition, (12. ip . This has been proved to all orders in the abelian case. For the non-abelian case 
the exponentiation for the Wilson loop for ordinary QCD has been demonstrated in [541 155] . 
Although we assume that that analysis extends to the Gribov-Zwanziger case we note that 
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for our one loop computation, the exponentiation does indeed occur which is sufficient for our 
present purpose. 

For the moment we follow the early approach of \36\ [371 EH] for the canonical formalism of 
QCD and the Lagrangian is given by 

= - ^G^G''^"' - ^{d^A^f - c^d^'D^c'' + iV^^^W*^ (2.5) 
where GJ^,^ is the field strength, the covariant derivative, D^, is defined by 

D,i;^' = + igT^jA'^^i;^' (2.6) 

g is the coupling constant and we have included the usual linear covariant gauge fixing term 
with parameter a. It is formally present to allow for the non-singular inversion of the quadratic 
part of the Lagrangian in momentum space to determine the propagators but we will set a = 
throughout thereafter. Associated with this are the Faddeev-Popov ghosts, c"" and c", and 
massless quarks, ip^^ , are included where the indices have the ranges 1 < a < Na, 1 < i < Nf 
and 1 < / < Np where Np and A'^ are the respective dimensions of the fundamental and adjoint 
representations and is the number of quark flavours. To recover the static set-up of the colour 
sources we take, [36 | [37 1 [38]. 

r^{x) = gv.T'^ [^(3) (x - ir) - J^^) _ ir')] (2.7) 

where = rj^Q is a unit vector and we will set r = |r — r'|. This vector will be present 
throughout our calculations and in effect projects out the time component of the gluon it couples 
to. Whilst we will be concentrating on the Gribov-Zwanziger case throughout, which relates 
to a theory of a confined gluon, we note that the group generator present in the source will 
be taken as being the fundamental representation initially. Though one can regard the heavy 
sources as being gluons in which case one would use the adjoint representation of the source 
interaction. The calculation of the one loop static potential is the same irrespective of whether 
the generator is in the fundamental or adjoint representations. Computationally differences will 
only arise in the two loop static potential, |43j . At this point it is worth briefly mentioning that 
in the Coulomb gauge, which Gribov also considered in \W\, it is the time component of the 
gluon propagator which receives attention in regard to conflnement. So the vector is in effect 
the bridge between the Landau gauge we use here and the Coulomb gauge results in relation to 
a confining potential through the gauge invariant Wilson loop. 

Given this formulation of the static potential, it is possible to evaluate V{r) either directly 
in coordinate space, (Ml 1171 [38], or in momentum space, [Ml [40l El [42]. For the latter V{r) 
emerges after a Fourier transform and in our conventions we take 

V{r) = J J^e^'^-'Vik) (2.8) 

where V(k) is the momentum space static potential. Performing the angular integrations since 
y(k) will only depend on the momentum length k then 

Vir) = ^2 f dk k'Vikf^ . (2.9) 

At this point it is perhaps apt to recall differing forms of V{k) and examine their implications 
for coordinate space. Although we will consider forms relevant for this article, a more general 
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library of mappings between spaces can be found, for instance, in [56j- Taking 



Vc{k) 



Ac 

fe2 



Vlik) 



A, 



2^2 



and Vh{k) 



Aj, 



{k 



2\3 



then 



Vc{r) 



Viir) 



Ai 



Ac_ 

Anr ' ' Sit 

respectively which result from the elementary integrals 



and Vb{r) 



dx ■ 



smx 



vr 
2 



dx ■ 



smx 



— and 
4 JO 



dx ■ 



Ab 3 

r 

967r 



sm j; 



vr 
96 



(2.10) 



(2.11) 



(2.12) 



respectively. So a linearly rising potential will emerge from a simple dipole term. However, it 
is amusing to note that an asymptotic linear behaviour as r — )• oo can also emerge from other 
forms of V{k) such as that given by including in Vi{k) any reasonable function of k, which is 
non-singular and does not upset the integral convergence, purely on dimensional grounds. 

In the early work of [36l EZl [38] the static potential was calculated at one loop directly in 
coordinate space. The presence of the static sources modifies the Feynman rules and introduces 
Heaviside step functions of the source interaction positions. Whilst it is relatively straightforward 
to compute the one loop potential, since we will ultimately be dealing with a formulation of QCD 
which is effectively massive, it is more natural to work in momentum space. The momentum 
space formalism was developed in [39l[40l|4Tl[42] and we defer to those articles for the justification 
of the more technical aspects rather than unnecessarily repeat them here. Suffice to say that the 
momentum space Feynman rules for the source gluon couplings do not alter from those recorded 
in an appendix of \41\ and we note that [41j gives comprehensive detail on many of the technical 
aspects of computing the static potentials. Moreover, the advantage of a momentum space 
approach is that it is in principle more straightforward to extend to two loops. For completeness 
in reviewing the earlier static potential result for pure QCD, we note that the one loop MS result 

of [Ml 133 EB] is 



^(p) 



IGTr^C^a 



QCD 



1 + 



31 11 



In 



Ca + 



In 



20 



TpNf 



(2.13) 



where fi is the mass scale introduced to ensure the coupling constant is dimensionless when 
dimensional regularization is used, a = ^(^/(IGTr^) and the colour group Casimirs are defined by 



Tr T^T 



ab 



rj-iClrj-iQ 



Cpl 



jacd jbcd 



ab 



(2.14) 



and / is the Np dimensional unit matrix. 



To extend the static potential to the Gribov-Zwanziger case we return to the definition (12. ip 
and ()2.3p . In Gribov's original work, [10], the key observation in the Landau gauge was that of 
the partition of configuration space into regions defined by zeroes of the Faddeev-Popov operator 
Ai(A) = — d^D^ which is hermitian. This is because one has copies of the gauge configuration 
satisfying the same gauge fixing condition. The region containing the origin, = 0, is referred 
to as the Gribov region and denoted by Gribov showed that to exclude copies from the gauge 
fixing procedure, the path integral measure must be restricted to the region which effectively 
introduced a cutoff into the theory. Consequently, the structure of the theory was altered in 
that the gluon propagator ceased to be of a fundamental form in the infrared region. This was 
primarily due to the appearance in the Lagrangian of an additional non-local term. Its origin is 
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in defining the boundary of fl by the no-pole or horizon condition and leads to the presence of 
a mass parameter, 7, known as the Gribov mass. This is not an independent object but must 
satisfy a gap equation derived from the horizon condition, [10]. One outcome of this is that the 
propagator of the Faddeev-Popov ghost takes a dipole form in the infrared. Whilst this is of 
the form for a linear confining potential in coordinate space, the Faddeev-Popov ghost is never 
directly exchanged between gluons or quarks. Given this position it is straightforward to see 
what one must do to determine the static potential in the Gribov picture. One simply returns 
to ()2.3p and restricts the integration measure to 



Z[J] 



VAf,VipVipVcVc exp 



(2.15) 



Implementing the horizon condition this equates to, |12] . 



Z[J] 



VA^VipV^VcVc exp 



Grib 



+ r^'Al 



(2.16) 



where 



L 



Grib ^ ^QCD _^ ^AT^^a^ 



AO- 



2g^ 



where d is the dimension of spacetime. 



(2.17) 



To proceed we now use the localization of the horizon non-locality introduced in [ 12\ [T3l 
[T71 [TBj by Zwanziger. This involves introducing an additional set of fields called localizing or 
Zwanziger ghosts, {0^^, and {uj'^,Ld'^}, where the former are commuting and the latter are 
anti-commuting. They are regarded as internal fields, similar to the Faddeev-Popov ghosts, in 
that they do not couple directly to quarks. The resulting path integral becomes 



Z[J] 
where 



/ 



VAfj, VipVipVc Vc Vp Voj Voj exp 



j d'x (lGZ 



+ r^'A'i 



(2.18) 



L 



GZ 



ab 



ucfp"'^ 



i^■iJabcJ^a^,^bc 



dNAl' 



and we have introduced the real and imaginary fields p'^^ and similar to |50] given by 



(2.19) 



^ ( ab , ■ t-ah 

7! + 



^ f ab 



ab 



(2.20) 



In choosing to work with the real fields /o^ and we will en route be correcting the error in 
earlier loop calculations, [32l[33], where the /o^^ propagator was erroneously treated in the initial 
computer algebra derivation. One main benefit is that (j2.19p is renormalizable, [18 ^ 130 1 [3T]. which 
allows us to do explicit calculations. We briefly note that in using (12.19P to perform computations 
the gauge symmetry is broken in much the same fashion as in the 7^ = situation via the (local) 
gauge fixing criterion in (j2.5p . However, the usefulness of ()2.19p in determining quantities of 
physical interest, in addition to renormalizability, resides essentially on two criteria. These are 
the gauge invariance of the object which is apparent, for instance, in the local or ultraviolet 
limit and the fact that 7 is constrained in the Gribov gap equation to be a (non-perturbative) 
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function of the coupling constant and is not an independent parameter of the theory. In other 
words (j2.19p has no meaning as a gauge theory unless 7 satisfies the gap equation defined by 
the horizon condition, [HI [12l [HI [HJ [EJ [HJ [13 [IHl [H]. For completeness we note the Landau 
gauge propagators of the gauge sector required for that computation and the current one are, 





-p)) = 


{Ai{p)eu\ 


-p)) = 


{Al{p)pl\ 


-p)) = 


{e^{p)c\ 


-p)) = 


{e^{p)p''X 


-p)) = 


{pf{p)pt'{- 


-p)) = 







■p,^Ap) 



where 



^ac^bd jabe jcde^i 



{u^\p)of{-p)) = - -^V,u (2.21) 

PAp) = V,u - ^ (2.22) 

is the transverse projector and we have a procedure similar to that discussed in [33J for handling 
the mixed propagator. The explicit renormalization properties of the Gribov mass and fields 
yO^*, ^"^^ and uj'^^ are predetermined by Slavnov- Taylor identities, [THl [301 [31], which have been 
verified by explicit computations in the MS scheme, [321 [33], in the ultraviolet situation where 
one can work with the 7 = limit. Those calculations were carried out with the symbolic 
manipulation machinery of FORM, [57], which we will also use here. The absence of a mixed 
P'^'^li propagator is due to the fact that the cross-term of the quadratic part of the Lagrangian 
is a total derivative which can be dropped, [SO]. However, we have retained the gluon /o^''-^^* 
vertex even though it can be written with a factor proportional to the Landau gauge condition, 
d^A^j^, [50]. Whilst this will vanish for transverse gluons in the Landau gauge we retain it as part 
of our Feynman rules since we will consider longitudinal corrections to the 2-point functions at 
one loop. 

Indeed given this, it is worth recalling how the propagators of (j2.2ip are constructed in 
practice since it is partly related to the situation with regard to other linear covariant gauges 
and justifies why we focus solely on the Landau gauge in the Gribov case. Also these comments 
are based around the more detailed analysis of [58], to which we refer the interested reader, 
where Gribov copies were considered for linear covariant gauges when the gauge parameter is 
small. First, we recall, [33], that in deriving ()2.2ip in order to avoid a singular determinant in 
Lorentz space in inverting the quadratic part of the momentum space Lagrangian, the gauge 
parameter is kept non-zero. The Landau gauge expressions, (I2.2ip . emerge when a is set to 
zero, [33]. Specifically, if we work in the basis of fields {^^, C^^, P^''} then with a non-zero a the 
matrix of quadratic terms in the Lagrangian in momentum space is 



^if'^'Hp) 



( -^2^"^ _ij2jacd Q 

_i^2fcab _p2§ac^bd g \ P^M 

\ _p2^ac^-bd 



+ 



/ _ 2_^ac _ ^^2 jacd 

_i"2fcab _p2§ac^bd g \ L^,^{p) (2.23) 

\ 
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where 

L,Ap) = ^ (2.24) 

is the longitudinal projector which satisfies the trivial relation P^u{p) + L^y{p) = rj^^. If we 
denote the matrix of propagators by uj^^'^^^p) then the propagators must satisfy 

f 6^P \ 

^{a6|cd}^^)n{cdbg}o-^^^) = | 6''P6'^'i | f]^^ (2.25) 

and the matrix on the right hand side is the unit matrix on the colour space for this sector of 
fields. Equipped with this the non-zero a propagators are 

{Alip)Aii-p)) = - [(^2)2 + C^y]^M-(P) - [(p2)2 + «C^^y]V(l') 

■ fabc 2 • fabc 2 

^ac^bd jabe jcde^A ^jabe jcde^A 

{e'{p)C\-p)) = - + p2[(p2)2 + C'A7l^'''^^^ ^ p2[(p2)2+^C^^4]V(P) 



{if{p)Pu{-p)) = 

{Pt{p)pf{-P)) = {u^^Hp)^fi-p)) = - . (2.26) 

These reduce to those of ()2.2ip in the a — )• limit. Though we have included the non-zero a 
forms since the longitudinal pieces play an important role in determining the one loop corrections 
to the gluon propagator. However, the propagators of (j2.26p have no meaning as such and must 
not be regarded as the propagators for the extension of the Gribov-Zwanziger Lagrangian to 
other linear covariant gauges. The Gribov-Zwanziger Lagrangian is purely a Landau gauge 
construction. To clarify this we recall certain properties in the construction of (j2.19p based on 
[58j . First, in the set of linear covariant gauges the Faddeev- Popov operator is only hermitian in 
the Landau gauge which allows for the classification of its eigenvalues into positive and negative 
values in that case only, [ini[IIl[l2l[ia[I3[l5l[ia[l6l[I71[l8l[l9]. This gives a natural way 
of dividing configuration space into definite regions. For other linear covariant gauges the loss 
of hermiticity means the eigenvalues can be complex and therefore there is not a natural or 
straightforward way of partitioning configuration space to even examine how Gribov copies are 
mapped between regions as one can do in the Landau gauge. Therefore, for non-Landau linear 
covariant gauges it does not seem clear how the path integral could be cut-off and the analogous 
first Gribov region defined by a no-pole condition as in the Landau case. Moreover, as the latter 
condition manifests itself as the non-locality to be localized in the Landau gauge one might 
simply assume that the natural extension to other linear covariant gauges is to merely use (j2.17p 
and proceed with a non-zero gauge parameter. However, this overlooks the fact that in the 
Landau gauge the gluon is transverse but in other linear covariant gauges it is not. Therefore, 
in the non-local term of (12.17P the gauge field in the denominator covariant derivative itself 
actually involves an additional non-local projection, |58j. Thus, for non-zero a the analogous 
Gribov-Zwanziger Lagrangian, (I2.19p . would have additional non-localities which would also 
need to be localized before any computations could proceed assuming any final local Lagrangian 
would actually be renormalizable, [58] . So, for instance, the intermediate propagators of (I2.26P 
that exist for non-zero a in deriving (j2.2ip have no true meaning. 
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3 2-point function corrections. 



Given that we are now using a version of the Gribov-Zwanziger Lagrangian which correctly 
involves the real and imaginary parts of the bosonic localizing ghosts, it is appropriate to discuss 
the one loop corrections to all the 2-point functions and propagators of the fields. Though we will 
concentrate primarily on the mixed sector for reasons which will become apparent later. Whilst 
the main details of this analysis parallels that given in [33], which was incorrect, we include it 
here for completeness. We define the 3x3 matrix of one loop 2-point functions formally by 

/ _p2^ac _^^2jacd q 





V 



+ 



+ 



+ 







.p2^ac^bd 




(3.1) 



with respect to the same basis {^^, where 



Qabcd 




Qabcd 


= Qp6''^6^'^ 


qL abed 


= Q^6''^6'"^ 


qL abed 





+ 
+ 



abe fcde 



ace rbde 



Wpr^f 



I jace jbde _j_ j^L/ jaoe 



+ Rd f 

L fabe rede 



7 



bde 



I j^L jabe jede 



represents the colour decomposition and 



jabed 
a A 



rp(brpc rpd) 



(3.2) 



(3.3) 



is totally symmetric and is the adjoint representation of the colour group generators. The 
quantities in the final two matrices of ()3.ip represent the one loop corrections and we have 



incorporated the fact that there is no one loop correction to the p^f-^^ 2-point functions. The 
longitudinal 2-point functions are indicated by the superscript L. Given ()3.ip we define a similar 
formal form for the inverse to one loop by 



jped 



17^ 



epq 



jedr jpqr 








[{p^y^+aCA-y"] 



jped 



V 



+ 





Q fcpq 
n^dpq 



l/ 



epq 



07 





jcdr jpqr 



y (jjped j^cdpq 



(Jjcpq 
J^cdpq 
j^edpq 



PtMu{p)a 



\ 







Lpv{p) 
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/ A^S'^P B^fP'i C^f'Pi \ 

^ fP^'^ E '^'^Pi Dp^Pi J 
for non-zero a, where 

j^cdpq _ J^^cp^dq _|_ jcpe jdqe _|_ q jcde pqe _|_ JJ^f^PI 



(3.4) 



j^cdpq 
jjalpq 
j^L cdpq 



D 



L cdpq 
j-^L cdpq 



D^S^PS'^'^ + J^fP^fdie + j^^fcdej^pqe ^ L^d^^PI 

DpS^Pb'^'i + JpfP^fl^ + Kpf^^fPI^ + Lpd'^^P'^ 

j^L^cp^dq _|_ jcpe jdge _|_ ^^L jcde jpqe _|_ JjL^cdpq 

J^L^cpgdq ^ jLjcpejdqe ^ jcde jpqe ^ ^L^f^q 

J^L^cpgdq ^ jLjcpejdqe ^ jcde jpge ^ LL^<^dpq 



(3.5) 



and the inverse satisfies (j2.25p again. The quantities in the final two matrices of (j3.4p represent 
the one loop corrections to the propagators. By multiplying out the elements of the left hand 
side of ()2.25p to one loop we can determine the relation of the one loop propagator corrections 
to those of the 2-point functions. For the transverse sector we find 



A 
B 
C 



1 



1 



{p'yx - 2iCAlVU - Ca7^ [Q^ + CaR^ + ^CaW^] 
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[{p^f + CAl^f 
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" [{p^y + Cai^] 
1 

, i'[Q^ + hCAW^] 
{p^?[{p^? + Cai^] 

, F = , G = 



nVx - {{p'Y - CAt)u + Z7V [Qi + CaRh + ^c^vrd 



{p 



,2\2 



P)2 



(p2)2 



p2[(p2)2 + CA7' 
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Rt) 





1 LP 



(3.6) 



to one loop. For the longitudinal sector, we retain for the moment the non-zero a and find 
formally, 



B^ 
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[{p^f + aCAl^f 



{p'fx'^ - 2iCAlVU'' - Cai'' + CaR^ + iC^Wf 



.,2 2ttL 



a 
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aV^ 
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(p2)2 



Li 



cL 
2^2 



(3.7) 



to one loop too. We note that the exact one loop expressions for each of these 2-point functions 
are recorded in appendices B and C. Whilst these represent our arbitrary a manipulations, we 
must restrict a to the Landau gauge for the propagators corrections to have any meaning in 
the Gribov-Zwanziger context. Since the one loop corrections are non-singular in a then taking 
a — >■ for the longitudinal sector we find 



Di 



(p2)2 
,2^2 



(p2) 



1^ 



2\2 
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,2^2 



[p 
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(p2)2 
(p2)2 



nL 



2\2 



oL 
^P 

2\2 



{P 



(3.8) 



Using the explicit values from appendices B and C one can in principle deduce the formal 
expressions for the one loop corrections to the propagators. For instance, given that = 
when a = then the one loop correction to the gluon propagator is actually transverse as is the 
mixed gluon ^'^^ propagator. Given the structures of D^"'^'^'^ and D^"''"^'^ then both the .^^^ and 
p^'' propagator corrections are not transverse like (|2.2ip . 

For the remaining fields and cjJJ, if we write the ghost form factors as 



then 



Dc{p^ 



{c\p)c\-p)) 



dUp^ 



p2 



{ojfip)Qt'i-p)) 



jac jfed -Pg^ (P 



-Vpu 



(3.9) 



1 + 



+ 



4p2 



tan 



p2 



8p2 

2 



+ 
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■ tan 
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8(p2)2 



1 + 



{p 



2\2 



Cao. 



1 + 



{p 



2\2 



+ 0{a^ 



Cai\ 
(3.10) 



Using the one loop Gribov gap equation, [10], then ^(p^) is 0((p2)2) p2 _^ q hence 



both propagators enhance. Likewise the propagator enhances whilst one colour component 
of the propagator also suggests there will be some sort of enhancement for that field too, 
j. We will return to this latter point in more detail in section 6. 



Equipped with these 2-point functions and the formal expressions for the one loop correc- 
tions to the {^Ji, C^^, pji^} matrix of propagators we can study the effective renormalization 
group invariant coupling constant. The behaviour of this effective coupling constant in the zero 
momentum limit has been the subject of debate over the years and whether it freezes to a finite 
non-zero or zero value is still unresolved. Defining the gluon propagator form factor as 



{Al{p)Al{-p)) 



p2 



-Ppuip) 



then the effective coupling constant is defined by 

a^fif(/) = asif^)DAip'){D,ip'))' 



(3.11) 



(3.12) 
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where we have introduced the more common strong couphng constant = /{A.tt). Using the 
zero momentum values for the 2-point functions recorded in appendix B we have 



«f(0) 



hm 



aip) l + CA(|ln 



Hence, evaluating this we find 



1 + 



^A(iln( 



«f(0) 



16 



+ 



8 ^ 8VCl7"(/^) 



o(/x) 



(3.13) 



(3.14) 



which is different to that of [33j since the effect of the omitted propagator has been correctly 
included. Numerically for SU{3) we have af{0) = 1.698 which is 4% lower than the value 
quoted in [33]. 



4 Static potential. 

We now return to the problem of computing the static potential. Given the reformulation of 
the basic QCD Lagrangian, ()2.5p . to incorporate the Gribov problem in terms of additional 
fields /9^^, u}'^^ and (D^*, it is important to stress that these localizing fields are regarded as 
internal and do not couple to the static potential source (|2.7p of (|2.4p . Ultimately one begins 
with a Lagrangian, (I2.19p . involving a gluon and it is this quantum which is confined. Thus the 
Feynman rules for the source coupling are the same as those for the usual perturbative case. The 
main difference is that one uses the Feynman rules for the Gribov-Zwanziger Lagrangian, (j2.19p . 
The latter is not a trivial point if one considers the Feynman diagrams contributing to the static 
potential. Concentrating, for the moment, on the single particle exchange graphs, the leading 
graph is illustrated in Figure 1 where the sources are represented by the thick vertical lines and 
the gluon by the spring. However, at the next order there are corrections to the exchanged 
particle and due to the mixing in the {^J^,^^''} sector, there are extra graphs aside from the 
gluon self-energy corrections. These are illustrated in Figure 2 where the central blob denotes 
all possible one particle irreducible one loop corrections. The split propagators involving gluons 
and denoted by the double line, are the propagators. One key point is that as there is 

no direct coupling of either bosonic localizing ghost to a source then there is no single exchange 
of this field akin to the diagram of Figure 1. Whilst the exact one loop corrections to all the 
2-point functions, which lead to the propagator corrections, are known explicitly, we still have 
to assemble all the pieces for the full static potential. For instance, in addition to the graphs 
of Figures 1 and 2, there are source gluon vertex corrections as well as double gluon exchange 
boxes at one loop. 



Ftnnnnnnnrdinr 



Figure 1: Tree contribution to static potential. 

Our calculation proceeds in general terms along the same lines as the earlier work of \39\ 
H0| H2] . The main difference is that Feynman integrals with source propagators as well as 
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Figure 2: Self-energy corrections to single gluon exchange graph. 



Gribov type propagators, as opposed to massless ones, have to be handled for non-single exchange 
Feynman graphs. For those where either the source propagator or vertex is corrected one has 
the potential problem of infrared divergences occurring. In the perturbative case the momentum 
space resolution of this was given in, for example, [321 113 E] • Briefly with massless propagators 
and a source propagator, l/{kv), an infrared divergence arises at the vertex correction but 
this cancels the double gluon exchange infrared divergence. However, this statement needs to 
be qualified since in dimensional regularization, which we use here, the regularizing parameter 
e, where d = 4 — 2e, cannot distinguish between infrared or ultraviolet infinities. Indeed in 
the original one loop Feynman gauge calculation of [37] they actually cancel in the one loop 
correction to each of the vertices of Figure 1. For other linear covariant gauges the infrared 
infinities cancel in the combination of double gluon exchange and source vertex corrections. In 
|41j this infrared problem was discussed in detail where a fictitious infrared regularizing mass 
was introduced to explicitly isolate the infrared divergences. They were then shown to cancel 
for an arbitrary linear covariant gauge. Indeed it was argued that at one loop these infrared 
divergences are expected to cancel since the full Wilson loop graph that the infrared divergent 
graphs each actually originate from, is the same diagram before it is cut open for the static limit. 
As that original Wilson loop topology is infrared safe the infrared cancellation could be regarded 
as a consistency check and the intermediate regularizing fictitious mass safely removed. In the 
Gribov case the presence of a natural mass arising from the Gribov parameter in the gluon 
propagators serves to act as a non-fictitious infrared regularization. Therefore, one can trace 
the infrared infinity naturally and verify its explicit cancellation. 

Having discussed this technical issue then all that remains is to describe the calculational 
set-up. The Feynman diagrams contributing to the one loop static potential are generated 
using the Qgraf package, ^59j . In total there are 1 tree graph and 31 one loop graphs which 
is more than the non-Gribov case due to the extra fields and mixed propagators. This latter 
total is broken down into one snail graph correcting the Gribov parameter of Figure 1, 18 
single exchange graphs, 10 source vertex corrections and 2 double (gluon) exchange graphs. The 
topologies which are additional to those of Figures 1 and 2 are provided in Figure 3. All the 
graphs in Qgraf notation are then converted into Form notation in order to harness the power 
necessary to handle all the resulting intermediate algebra. This conversion adds all the colour 
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and Lorentz indices to the fields. The Feynman rules are automatically inserted to produce the 
set of Feynman integrals which need to be evaluated explicitly. For the single particle exchange 
diagrams this amounts to applying the basic integrals which were already treated in [33] where 
the exact one loop corrections to all the 2-point functions in the theory were given. 



Figure 3: Additional one loop topologies for static potential. 



The final part is to handle the remaining box and source vertex corrections which will involve 
source propagators. Therefore, we define the general scalar master integral 



Ii{a,P,X,ml,ml,p,v) 



d<^k 



1 



(4.1) 



for arbitrary masses mi and m2 where is the external or exchange momentum and will 
satisfy pv = 0, with v"^ = 1, in the static limit throughout. Such integrals emerge after repeated 
substitution of the elementary relations 



kp 



k^ + p2 



{k-pf 



and 



1 



1 



1 



(4.2) 



(4.3) 



Although (j4.ip has been given in [41j . for completeness we give several intermediate results we 
used, in our notation. Using integration by parts it is possible to write a recurrence relation for 
(liTl) as 



All (a, /3, A + 1, m^, m2,p, v) 



— 2ali{a + 1, /?, A — 1, m\^m\^p, v) 

— 2/3/i(a, /3 + 1, A — 1, m\,m\,p, v) 



(4.4) 



for A > 2, which allows one to reduce integrals with two source propagators to a standard source 
free integral which has already been considered in the Gribov case in [33J. For an odd number 
of source propagators, this reduction will eventually lead to one source propagator but it can 
be treated using elementary contour integration in the time integral. So, for instance, we have. 



h{W,mlmlP.v) 



(27r)^-i [k2 + m2][(k - p)2 + ml\ 



(4.5) 
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In other words in this case the treatment of the source propagator reduces the dimensionahty of 
the original integral by one. Finally, any remaining source propagator free integral with prop- 
agator powers greater than unity can be reduced by standard integration by parts identities. 
One final aspect of the automatic calculation concerns the internal renormalization of the pa- 
rameters and sources. We follow the standard procedure of [60j for this and compute the static 
potential initially as a function of bare parameters. The renormalized variables, and hence the 
underlying counterterms, are introduced by replacing the bare parameters by the renormalized 
ones with the explicit Landau gauge renormalization constants included. This also includes the 
renormalization associated with the source itself and that renormalization constant is derived 
in the same way as discussed in [3^ but for the Landau gauge. For completeness, we note that 
the relevant MS Landau gauge renormalization constants, in our conventions, are 
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1 + 0{a^ 



(4.6) 



where Jq^ = y/Zj J"'^ relates the bare and renormalized sources. Obviously a non-trivial check 
on the renormalization procedure and the implementation of the explicit values for the master 
Feynman integrals is that the ultimate static potential is finite and no poles in e remain when 
all the contributing diagrams are assembled. We note now that this is indeed the case. 

Having described the underlying tools for our calculation we now record that the one loop 
static potential for the Gribov-Zwanziger Lagrangian in the Landau gauge is 
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where we use the notation that in four dimensions Pfj, = {po,p) and we have introduced the 
intermediate functions ?7i(p^) and r]2{p') for compactness where 
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and 
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There is one main check on this result aside from the finiteness one above. This is that we 
recover the usual perturbative result in the t^ ^ limit. In other words we find 
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where we note that ()4.10p agrees exactly with the earher result of [Ml EZl . Whilst that was 
originally computed for the Feynman gauge taking the limit where the Gribov mass is removed 
recovers the perturbative result in the Landau gauge and its agreement with a gauge independent 
result is another non-trivial check on our computational set-up. 

We can now examine (j4.7p in the zero momentum limit by expanding the exact result in 
powers of and find 
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where the order symbols denote the higher order terms in the momentum expansion and two 
loop corrections separately. Thus (j4.1ip would imply that 
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as a first approximation with no emergence of a dipole. To try and understand the implications 
of this freezing for the one loop coordinate space potential we first consider the leading term 
of (j4.1ip . For 7 = the exchange of Figure 1 leads to the usual Coulomb potential but for 
non-zero 7 it is not clear what the _). g limit of the leading term of ()4.1ip implies for the 
corresponding r — t- 00 limit. If we write 



(4.13) 



n=0 



where the subscript n labels the loop order, then the first term of (j4.7p gives 



^o(p) = - 
whence the Fourier transform, (j2.8p . leads to 
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The presence of the exponential factor is not unexpected since it is reminiscent of the Yukawa 
potential for a massive field. However, the Gribov situation effectively corresponds to a width 
which is reflected in the trignometric factor. Now examining the limit corresponding to — >• 0, 
which is r — )• cxD for real r, we see that Vo(r) — )• 0. This is consistent on dimensional grounds 
with naively taking the Fourier transform of the leading term of (j4.1ip . As an aside we remark 
that similar potentials of this Friedel form have recently emerged in leading order in models of 
plasmas and stellar nuclear reactions, |6 H 162 1 [63 l 164 1 165j . Briefly the underlying theory is based 
on a Higgs-like effective Lagrangian. By considering perturbations about a background then the 
gauge field fluctuations develop a Gribov type propagator. One claim is that a linearly rising 
potential will emerge in a particular limit if the effective Higgs mass term has the wrong sign. 

Turning to the one loop correction we do not have the luxury of being able to take the 
full Fourier transform of (j4.7p analytically. However, by the same heuristic argument as leading 
order, the freezing to a finite value, (|4.12p . leads to a similar large distance limit for our potential. 
In other words V{r) — )• at one loop and there is no linear growth. However, on dimensional 
grounds such a term will clearly not be Coulombic in this limit. A 1/p^ term would be required 
for that in the zero momentum limit. Indeed given this the next step would be to try and 
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understand how the Liischer term, [66], emerges in the Gribov-Zwanziger formahsm as r — )• oo. 
Within the present context it is not clear how such a term could appear. For instance, aside 
from the fact that the analysis of [66] producing this term was based on chromoelectric flux 
tubes, the coefficient of the Liischer term is universal and independent of the coupling constant. 
One prospective way this could happen is that the coefficient of a 1/p^ term, when it arises, is 
fixed by the gap equation satisfied by 7 with or without some renormalization group running of 
parameters to some non-perturbative fixed point. Moreover, the Liischer term appears together 
with a linearly rising potential term, [66], which would originate in momentum space from 
a dipole. However, a dipole behaviour in the zero momentum limit is clearly not evident in 
()4.7p . Although we have not taken the complete Fourier transform of (14. TP . we have examined 
it in detail for a confining behaviour to understand why a dipole is absent and there are some 
promising features. 

Clearly in (14. 7p there is an explicit l/(p^)^ term which gives a stronger singularity than that 
of the usual Coulomb case. Indeed concentrating on this term for the moment and setting 

' 61447r ^ ' 



in (I2.10p we would have 



Vi{r) = - ^^^^ ^/ r (4.17) 
^ 491527r2 ^ ' 

which would actually be a linearly decreasing potential. Ignoring this point for the moment, 
there is, however, another dipole singularity in (j4.7p and this occurs in such a way that there is 
no overall singularity as — )• 0. The specific term is 



30727r2(p2)2'^'^ 



p2 



(4.18) 



and in the zero momentum limit the inverse tan function tends to tt/2 to cancel the previous 
pure dipole we considered. Thus, a linearly rising potential will not emerge. Actually for (|4.18p 
it is in fact possible to take the Fourier transform (12. 8p and verify that the overall absolute 
asymptotic behaviour as r — t- is linear and precisely cancels that from the pure dipole term 
Fourier transform. However, we recall what was noted in [33] concerning the overall sign of 
7^. In studying the zero momentum limits of (|4.16p and (j4.18p we have tacitly assumed that 
7^ is positive. The sign of 7^, though, is not predetermined in (I2.19p . It can only strictly be 
determined in relation to other evidence. Indeed it was suggested in [33] that for consistency 
with other results, such as the power correction structure of the strong coupling constant in the 
next to high energy limit, that 7^ was negative and our sign convention needed to be modified 
by mapping 7^ — — 7^ in situations where results depended explicitly on 7^. It is worth noting 
that this is not cause for alarm since a similar situation always occurs with the choice of the 
sign of the coupling constant, g, in the covariant derivative. Its sign cannot be determined since 
in any computations it appears in the combination g^. In other words the choice of sign of g is 
a convention as is that of 7^ with the difference being here that one can actually make contact 
with other methods, such as the strong coupling constant power corrections, to determine it, 
[33]. In our current set-up we had not made any a priori choice of sign and so there is a degree 
of freedom to make a specific choice. Therefore, if we return to (I4.16P and ()4.18p . then since the 
former is an odd function and the latter is an even function of 7^ , flipping the sign of 7^ would 
mean that in fact the terms add and there is a net dipole. This would lead to a linearly rising 
potential. Although this is encouraging and would very much be consistent with the expectation 
that the Gribov-Zwanziger Lagrangian describes a confined gluon since, for instance, it satisfies 
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the Kugo-Ojima criterion, [MIISS]) it would on the other hand destroy gluon suppression, [16]. 
This can be seen if one examines the exphcit one loop corrections to the gluon propagator given 
in appendix B. More specifically examining the relevant pieces there, the net dipole arises purely 
from the 2-point function corrections in the one loop gluon propagator. Indeed this would 
be consistent with the Zwanziger localizing fields dominating the infrared as one moves towards 
the Gribov horizon. If such a scenario is to be credible, however, one needs to check that at 
next loop order a triple pole in does not arise in the static potential. Such a term is possible 
on dimensional grounds, for instance. We stress, though, at this point that these are merely 
interim speculative remarks on the consequences of the sign choice of 7^, which would require 
deeper consideration, and we will omit reference to this point hereafter. 

Therefore, returning to our original 7^ convention this dipole cancellation is not unexpected. 
Indeed if one examines the explicit exact one loop corrections to each of the transverse parts 
of the 2-point functions in the {^^, ^JJ^, p^^} sector, this dipole structure is present in each case 
in individual terms but with no net dipole. Thus the apparent dipole behaviour is not truly 
present overall. However, on more general terms this lack of a linearly rising potential should not 
emerge in the present set-up. Although we have followed an inherently perturbative approach the 
main lesson from Gribov's original article, which also underlies Zwanziger's construction, is that 
infrared properties emerge when the gap equation is realised. Then one is dealing with a gauge 
theory. For example, the gap equation leads to Faddeev-Popov ghost enhancement and more 
recently produced a qualitative non-zero value of a renormalization group invariant coupling 
constant in the zero momentum limit, (13.14^ . However, at no point have we implemented the 
gap equation within our static potential. (The simple replacement of 7 by a non-perturbative 
function of a is not what we mean by this here.) Somehow, one would expect that a linearly rising 
potential can only emerge when the gap equation is used explicitly. In a later section we will 
present some considerations towards this point in the context of the static potential considered 
here. Though, recalling an earlier remark, at higher loop orders there would appear to be 
nothing in principle preventing pure terms such as l/(p^)'^ and higher emerging on dimensional 
grounds. These would have to have a compensating piece to exclude power terms in r remaining 
in the coordinate space potential. Of course such terms would be irrelevant if the sign associated 
with it indicated a decreasing contribution to the potential in contrast to an increasing part. 
Alternatively such terms might resum in such a way that there is matching to a linear type 
potential. One final point concerning the linear term relates to the soft BRST breaking of (I2.19P 
arising from the 7^ term, [H [HI [H |30] . There is an understanding that for gauge independent 
quantities results should depend on 7^^ and not on 7^. However, on dimensional grounds a linear 
term in a potential must be accompanied by a factor of 7^ and therefore it seems impossible 
to avoid having a 7^ dependence in the ultimate static potential which is derived from the 
underlying gauge independent Wilson loop. Though, of course, any such scale can eventually in 
principle be related back to A^^. 

Next we make some specific remarks concerning (j4.7p . Given that the gluon is now not 
a massless entity in (I2.19p . we note that (14. 7p has a pole and physical cut at = 2\JCa')^ ■ 
This is a threshold type feature due to massless fields with a non-zero width. The associated 
denominator factor in the term with this cut arises from the Gram determinant of the integration 
by parts reduction formula for the one loop Feynman integrals without source propagators. In 
[13j . the glueball correlation function was analysed as a spectral density function and it was 
argued that there was a bound state of the same mass squared value, 2\JCa1^ ■, ^-^id it was 
identified as a potential glueball state. It is amusing that there is a threshold effect at the 
same mass value in the static potential which is in some sense a scattering amplitude. It is 
not clear, though, if this corresponds to a stable bound state. However, it is worth recalling the 
computation of |67j where the glueball spectrum was determined by including a massive gluon in 
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the evaluation of the leading order gluon scattering amplitude. This was then used to construct 
a potential within which bound states could be formed. However, the scattering amplitude with 
the massive gluons did not contain an explicit confining potential but this was added by hand 
prior to solving for the spectrum. Though strictly speaking the additional confining piece not 
only incorporated a linear rising part but also modelled string breaking. This final form of the 
potential energy appears to have endowed a stability on the glueball states. Therefore, it does 
not seem unreasonable to expect that bound state stability requires an explicit dipole term in 
(j4.7p . with string breaking, for similar reasons. For completeness we note that in [3T] massive 
W and Z gauge bosons were included in the static potential formalism for a standard model 
study but there the real mass did not lead to the same threshold singularity as ()4.7p . Next we 
note that one can isolate the contribution to ()4.7p from purely single particle exchange diagrams 
corresponding to all the graphs of Figures 1 and 2, since it too has similar dipole behaviour. 
Indeed it could be the case that the single exchange graphs have a net dipole term which is 
cancelled by a matching piece from the box and vertex graphs. However, analysing the single 
exchange contribution in the same way as ()4.7p the zero momentum limit is also non-singular. 
Though in comparison with (14. 7p there is an additional dipole-like term. As — t- these three 
terms conspire in a similar way as before to exclude an overall singularity. 

Although the dipole term does not emerge, we close this section by discussing other aspects 
of (|4.7p in the context of [36l |37l [Ml EH [40l |4l]. In the perturbative approach one can define 
a strong coupling constant which matches the usual ultraviolet behaviour of a strong coupling 
constant in the large momentum limit. It is based on the momentum space behaviour of the 
potential and is defined by 

^(P) = - ^«v(p) (4.19) 

where the subscript y denotes the F-scheme, [681 EH]. One advantage of this definition is that it 
is derived from a gauge independent quantity which is effectively the force between two coloured 
objects. The low momentum behaviour can be studied. Clearly in the Gribov-Zwanziger case 
at both leading order and one loop ay(0) = 0. Though for the former the behaviour in the 
zero momentum limit is O ((p^)^) but 0{p'^) in the latter. If, however, an overall dipole piece 
dominated the infrared then ay(p) would be singular at low momentum when one factor of 
is cancelled from the coupling constant definition. This would be consistent with confinement 
in this scheme when the strength of a suitably defined coupling constant increases to produce 
infrared slavery. Though in the absence of quarks there would appear to be no string breaking 
related to saturation discussed, for example, in [67]. Given these comments it might be an 
interesting exercise for lattice or Dyson Schwinger methods to be applied to the static potential 
specifically in the Gribov-Zwanziger Landau gauge case. This may be particularly apt given the 
current debate over whether the scaling or decoupling solution is the correct picture. Briefly, in 
essence we have concentrated on using the formalism associated with the scaling solution. In the 
decoupling scenario, [22l[23l[2H[25l[26l[271[28l[29j, it transpires that the gluon is not suppressed in 
the infrared, since the propagator freezes to a finite non-zero value but the form factor vanishes, 
and the Faddeev-Popov ghost is not enhanced. It has a behaviour which is close to a 
dependence as ^ rather than the dipole of the Gribov scenario, [22l[23l [2il [25| [26l[27l[28t 125] . 



5 Gap equation. 

In this section we focus on several aspects of the Gribov gap equation which defines the Gribov 
parameter 7 using ()2.19p to correct the error in the earlier two loop gap equation of ^32j. In our 
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current notation the gap equation is defined by the expectation value 



(5.1) 



Whilst the new features we present here are not immediately related to the static poten- 
tial, we believe they are important and the computations do contribute to understanding fur- 
ther underlying properties of the Gribov-Zwanziger Lagrangian. Throughout our calculations 
7 appears explicitly but is not an independent parameter of the theory. As emphasised in 
[m HI im [HI [151 [13 [13 [S]) the Lagrangian can only be interpreted as a gauge theory when 
the gap equation is satisfied explicitly. The original one loop expression was derived in [10] in 
the MS scheme and it is straightforward to invert it to determine 7 as an explicit function of 
the coupling constant. For instance, we have, [TO] . 
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which gives 
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exp 



or expressing it in numerical form we have 
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a + 0(0^) 
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(5.2) 



(5.3) 



(5.4) 



(In this section we revert to the MS scheme rather than continue with the ^/-scheme.) This 
clearly demonstrates the non-perturbative behaviour of the Gribov mass parameter. However, 
as the two loop MS correction to the mass gap is now available, [32], one interesting question 
is whether it is possible to rearrange the explicit expression to preserve the non-perturbative 
property emerging from the one loop result. Therefore, to do this we take 



co[l + ciCAas{n)]exp 



CAOisi^j) 



(5.5) 



as an ansatz for a corrected form of (j5.3p . Substituting this into the explicit gap equation, whose 
correct form now is, 
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where S2 = (2-v/3/9)Cl2(27r/3) and Cl2(x) is the Clausen function, it ought to be possible to 
determine explicit expressions for the unknown coefficients 60 and Cj. If one cannot find a 
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solution then an alternative approach would be needed. Remarkably, it transpires that to two 
loop order one can achieve the inversion for arbitrary numbers of massless quarks and we find 
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To assist with comparisons we have evaluated the coefficients numerically for both SU{2) and 
SU{3) colour groups and for several values of and these are given respectively in Tables 1 
and 2. Finally, we note that to this order in the gap equation ()5.5p can be rewritten as 
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Co exp 
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(5.8) 



with the same values for cq, ci and bo. 



Nf 


^0 


Co 


Cl 





16.913 


9.052 


- 0.133 


2 


18.383 


8.994 


- 0.180 


3 


19.303 


8.963 


- 0.210 



Table 1. Numerical values of parameters for SU{2). 



Nf 


bo 


Co 


Cl 





16.913 


9.052 


- 0.133 


2 


17.864 


9.014 


- 0.163 


3 


18.383 


8.994 


- 0.180 



Table 2. Numerical values of parameters for SU{3). 



Although the perturbative correction to the one loop form appears to be converging the 
actual convergence cannot fully be commented on until a three loop gap equation is available. 
One indication of the issues related to this is the change in value of cq between loop orders. In 
one sense this is not unexpected since the ansatz we have taken has a relatively simple form 
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and may not be the best approximation to the full relation for 7 as an explicit function of the 
coupling constant which is not known. To illustrate the subtleties of solving the gap equation 
at two loops, we had tried to solve ()5.6p numerically as a quadratic in the coupling constant. 
However, this was problematic since for certain values of the parameters, it resulted in solutions 
where there was an imaginary part. This was in contradiction with the understanding that 7 
and the coupling constant are real parameters. Moreover, if the solution to the quadratic had 
turned out to be real, it was not clear whether this would persist at higher loop orders when one 
would in principle have to solve cubic and quartic equations. The potential proliferation of roots 
would at some point be problematic and render such an approach meaningless. Therefore, we 
took the ansatz approach discussed above. As an aside we note that with the new two loop MS 
gap equation, (j5.6p . we have verified that the Kugo-Ojima confinement criterion, [Ml [35], still 
holds and that both the Faddeev-Popov and w^'' propagators enhance in the zero momentum 
limit at two loops, as before. 

We turn now to another aspect of the gap equation and that is to do with its form in another 
renormalization scheme. Whilst the MS scheme is standard in (high order) loop calculations 
it is not always the most appropriate since it is not a scheme motivated from a physical point 
of view. One such scheme, however, is the MOM scheme, [70^ I71j . which has also been used, 
for example, in lattice calculations of the gluon 2-point function and the evaluation of the 
renormalization group invariant effective strong coupling constant, |72[ [73| I74j . The latter is 
derived from the gluon ghost vertex, though the triple gluon vertex has also been considered. As 
the MOM scheme seems to be the one most studied in that context it would seem appropriate to 
develop the Gribov gap equation in the MOM scheme to see if there is any scheme dependence. 
Therefore, we repeat the computation of section 3 but for the MOM scheme which is not a 
trivial exercise. 

As a starting point we note that the exact expressions for the one loop 2-point functions are 
required as a function of 7. These are necessary since the MOM scheme is a mass dependent 
scheme and therefore the full expressions are needed to determine the correct finite parts for the 
MOM renormalization constants. However, as in the 7 = case the MOM renormalization has 
to be carried out in such a way that the underlying Slavnov- Taylor identities are not violated. 
For the Gribov-Zwanziger Lagrangian specifically, this means the renormalization constants for 
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Z^ and Z^, must satisfy, [TSl [30l [3T|. 
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where Zg is the coupling constant renormalization constant. Given this we have proceeded as 



follows. First, we determined the wave function renormalization constants for A'^, and by 
rendering the respective 2-point functions exactly equal to unity when the external momentum 
satisfies = fj? where ^ is the mass scale of the renormalization point. We record the explicit 
forms for Z^, Zp and Z^ in the MOM scheme are 
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and 
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We have checked that the latter is consistent by renormalizing the Faddeev-Popov ghost 2-point 
function in the same way and verified explicitly that the first of the Slavnov- Taylor identities 
holds at one loop. Equipped with these renormalization constants then that for 7 is now already 
fixed in the MOM scheme through the second identity of (j5.9p . In order to verify that this is 
not inconsistent we have repeated the one loop gap equation computation but using these MOM 
renormalization constants. This is achieved in the same way as ^32j, by closing the legs on the 
mixed propagator which gives the one loop contribution to the gap equation. As this is a vacuum 
diagram there is no external momentum to set to a mass shell value to extract the MOM scheme 
expression for which is why we have proceeded with the 2-point functions first. Therefore, 
the renormalization constants which have already been determined will remove the divergences 
and their finite parts will infiuence the final form of the gap equation. Ultimately we find the 
relatively simple expressiorQ 
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where for later purposes we have introduced a shorthand definition of the right hand side of the 
gap equation. Given that we are following the standard way to proceed in the MOM scheme there 
could be a doubt as to whether this is ultimately correct. However, in the Gribov-Zwanziger 
context the main check on this is whether the Kugo-Ojima criterion of Faddeev-Popov ghost 
enhancement, [34^ [35], emerges in the c° 2-point function in the zero momentum limit when 
the MOM gap equation of (|5.12p is satisfied. We note that not unexpectedly it does so. As a 
final point we have returned to the renormalization group invariant effective coupling constant 
defined from the gluon ghost vertex and reconsidered it in the infrared limit similar to section 



*This corrects a sign error in tlie expression given in [75j 
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3. Using the MOM versions of the respective form factors, D^'^^{p'^) and D^^^{p'^) then we 
define the MOM effective renormahzation group invariant couphng constant as 



MOM 



(5.13) 



Equipped with this we find that the couphng constant freezing calculation becomes 
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where we have displayed the rationals in the numerator explicitly in order to compare with the 
analogous MS computation of [33]. Hence, 
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16 



(5.15) 



This is the same value as the MS scheme, [33J, as expected and can be regarded as a robust check 
on the finite parts of the MOM renormahzation constants we derived. One final comment on 
the MOM gap equation in the context of the earlier discussion is that it is evident that inverting 
()5.12p is not as straightforward as that for the MS scheme and we have not proceeded. This 
is partly because the two loop expression is not available in order to do a proper comparison 
but also because if it was, it is clear that nothing substantial would emerge from a complicated 
explicit expression. Indeed a two loop calculation of the MOM gap equation would require the 
finite part of the 2-point functions exactly which would be a huge computation. To elaborate on 
the complexity of this problem in order to obtain, for instance, the wave function renormahzation 
constants there are well over 1000 two loop Feynman diagrams correcting all the A'^, p'^^ and 
^Jj'' 2-point functions. Aside from this, it is not clear if all the master scalar two loop self-energy 
Feynman integrals have been evaluated for all the necessary massive propagator configurations. 
For instance, the exact expression for the scalar master two loop self-energy corrections with 
purely 3-point vertices with all masses equal is not known. Even it were available one would still 
require the expression for different masses in order to determine the expressions for the Gribov 
mass combinations of ^i^JCAl"^- Indeed to appreciate how complicated this problem would be, 
one can examine the structure for the simpler two loop sunset self-energy correction for arbitrary 
masses, discussed in [76j . The arbitrary mass expression in this simplest case involves Lauricella 
functions where the arguments are a function of the masses and the external momentum. In 
an MOM renormahzation such functions, and the as yet undetermined master integrals for the 
two loop topologies with four and five propagators, would have to be evaluated exactly at the 
point where p^ = p?. Therefore, given these considerations it is inconceivable that a two loop 
MOM gap equation is viable in the immediate future, even if such a quantity were of interest. 
Finally, we note that the different values for the renormahzation group invariant and ^-scheme 
coupling constants at zero momentum are not inconsistent. In this region there is no unique 
way of defining the strong coupling constant. 



6 and enhancement. 

Whilst our potential has been derived from the usual static potential formalism at one loop, it 
is clear that it lacks a linearly rising part. This is due to the cancellation of similar l/{p'^)'^ type 
terms as p^ — )• 0. However, from experience of the Gribov-Zwanziger Lagrangian any infrared 
behaviour which is in qualitative agreement with expectations from other work has always 
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required the explicit use of the gap equation satisfied by 7. For example, the enhancement of 
the Faddeev- Popov ghost is due to the gap equation, [10], and we recall that the theory can only 
be regarded as a gauge theory when the gap equation is satisfied. Indeed the ghost enhancement 
is also a key component of the Kugo-Ojima confinement criterion, |341I35). in Yang-Mills theory. 
Its role in the Gribov-Zwanziger path integral has recently been clarified in [77j and is also used 
as a boundary condition in Dyson Schwinger studies, [29]. Whilst the Faddeev-Popov ghosts 
and localizing ghosts ujI^ satisfy the enhancement, they cannot play a direct role in actually 
confining a gluon since they are Grassmann fields. Hence they cannot be directly exchanged in 
gluon-gluon interactions as is evident from the graphs of Figures 1 and 2. Instead one clearly 
requires a commuting field to enhance to be at least in a situation where a singularity 
could be exchanged. Unlike the original expectations of Mandelstam and others, [3 [3 |8l [9] , 
this does not appear to be the gluon in the Gribov-Zwanziger set-up. However, as recently 
pointed out by Zwanziger there appears to be a clue in the enhancement of the propagators of 
the bosonic ghosts from a Schwinger Dyson analysis, [50]. Therefore, our aim in this section is 
to first demonstrate that that enhancement can also actually be accessed in perturbation theory 
and then discuss its implications for the static potential. Aside from the fields c° and uj'^ the 
explicit forms of the 2-point functions for and have the potential for enhancement. More 
concretely the explicit one loop correction to the colour channel of the Lagrangian kinetic term 
for each of ^^'^ and p'^ , given in appendix B, are exactly equivalent to that of the Faddeev-Popov 
ghost 2-point function. To illustrate this for completeness here we recall the — t- behaviour 
of the 2-point functions are, 
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where to avoid confusion with the propagator we have formally indicated the inverse. The first 
terms on the right hand side of ()6.ip and (16. 2p clearly correspond to the one loop gap equation, 
(j5.2p . which would imply enhancement similar to the Grassmann ghost fields and the emergence 



of an infrared dipole form. The inversion of the p'^ 2-point function to deduce the propagator 



is similar to that of the Faddeev-Popov ghost and and so p'^ will have an enhanced infrared 
propagator too. However, a enhancement is not as straightforward to observe as those fields 



due to the extra terms in (j6.ip and the group structure. 



To examine how enhancement emerges, we first recall the situation in the Faddeev-Popov 
case. There one first computes the 2-point function in the zero momentum limit and applies 
()5.2p . This produces a leading term of O ((p^)^ 



which one then inverts to discover the dipole 

ab 



infrared behaviour of the Faddeev-Popov ghost. Turning to the case the algorithm is the 
same but not as straightforward due to the mixing in the quadratic part of the {^uj'^u^} sector 
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of the Lagrangian. So the analogous inversion has to involve the full 2x2 mixing matrix. In 
section 3 this was carried out formally at one loop but without first enforcing the gap equation. 
Given the potential for enhancement of j50j in a Dyson Schwinger analysis, we reconsider 
the formal matrix inversion in a more general way. Specifically we will focus on the 
sector and in particular the transverse piece. We do this because at one loop V = and so 
the transverse part of ()3.ip becomes block diagonal. It is the upper 2x2 matrix which is of 
interest. If we now define this 2x2 matrix by Ag''^'^'^^, where we can drop the Lorentz indices, 
then formally 

where we use the more general decomposition 

Qabcd ^ Q^S'^c^bd ^ y^^jacejbde ^ ^^jabe jcde j_ g^^abcd j_ p^§abgcd ^ T^^a^^bc (6.4) 

and the quantities X,U, Q^, W^, TZ^ and 5^ are the formal 2-point functions including the part 
from the quadratic part of the Lagrangian. We have used similar notation to sections 2 and 3 
but in calligraphic font to indicate that these are not solely the one loop corrections. We have 
also included two extra terms, and 7^, to complete the basis. Although these are zero up 

to and including one loop they are required here since we will be multiplying A^'^''''^'^^ by its 
full inverse rather than drop the O(a^) part as was carried out in deriving (13. 6p . Here there 
will be extra group structures when, for example, two d'^'^'^ tensors are partially contracted. 
Whilst we are aiming at being general our matrix can only really be regarded as one loop in one 
sense since the colour group structure in the final element may not be the most general. Our 
choice there is motivated by what actually emerges from the explicit one loop computations. 
For instance, if quarks are present in higher loop diagrams then the tensor d'p'^'^ could also be 
present due to light-by-light subgraphs or some peculiar contraction of this with other tensors. 
For a comprehensive discussion of potential high rank (adjoint) tensors, see |78j . At present we 
make no assumptions about the behaviour of ()6.3p in the p-^ — )• limit but note that like Qp in 
(16. 2p . will be the key function in driving any enhancement. Also the loop order will play a 
key role in the inversion and we note that 



X = U = = 0(1) , = = = 0{a) , = = 0{a^) . (6.5) 

3n Ag"'*''^'^^ we define ' 
the same formal way by 



Given we define the general inverse Il2^'^^^'^\ which will be the matrix of propagators, in 



(6.6) 



where now 



^cdpq ^ jy^^cp^dq ^ j^jcpejdqe ^ J^^fcdejpqe ^ ^^^cdpq ^ ^^gcd^pq ^ (g^) 

similar to (13. 4p but allowing for the extra colour tensors to have a basis. For clarity we note 
that the two matrices must satisfy the standard inversion on the smaller subspace given by 

Ar'^'nl'"'"' = ( T ] (6.8) 



where the right hand side is effectively the unit matrix. Multiplying out the matrices explicitly 
leads to the formal linear equations satisfied by the 2-point functions and propagators. We have 

1 = AX + CaUB , Q = XB + (v^-Mi. + Ca1C^ + \caJ^U 
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QiPf: + h2C^W^ + h2Si:Ji. + a2S^C^ + T^M^ 



+ CaW^ + -C'aSi, + NaV^ + T^j + [CaJ^ + -CiC^ + V^+ M^j 
b2C^W^ + h2S^J^ + a2S^C^ + Q^M^ + T^P^ 

UB + Q^IC^ + ^CaW^J^ + ^CaW^/C^ + ^CA7^5J5 + CaTZ^JC^ - bsW^C^ 



(6.9) 



where the coefficients Oj and 6j derive from the group decompositions defined and discussed in 
appendix A where their exphcit forms are given for an arbitrary colour group. It is clear we have 
nine equations for nine unknowns. So it is a straightforward exercise to determine the general 
solution. First, we record that 



A 



[Qg + CaTZ^ + ^g^Wg] 
[(Qg + CA-JZi: + ^CaW^)X - CaU^] 



B 



U 



[(Qg + CaH^ + \CaW^)X - CaU^] 

(6.10) 

where we have assumed 7^^ = 7^ = initially in accordance to what we found at one loop. The 
explicit forms of the form factors for propagator are cumbersome. So for these cases we 
record the SU (3) expressions where the values for Oj and bi of appendix A have been used. We 
have 



2Qg 



^ ^3{3S^ - 2Wg) (5^ + 2W5) (5^ + W^) + 8(75^ + 3^^)21 
+ 16Q| + 2(275| + 205^^5-8^1)25 

[2Q5 + 355 + 3W^]-^ [2Q5 + 355 - 2W^r^ [2Q5 +S^ + 2W^r^ 
AW, 



[2Q5 + 355 + 3W^][2Qi: + 3S^ - 2W5] 



1 

2^ 



(4(35^ - W^)Q^ + 3(35^ - 2W^){S^ + W^)){'2U'^ - W^X - 2n^X) 

2Q^X + m^X -6U^ + 3W^X'\ [2Q5 + 355 + 3^5]"^ [2Q5 + 355 - 2W^]~^ 
4 [2Q^S^ + (355 + 2W^)iS^ - W5)] 

X [2Q5 + 355 + 3^5]"^ [2Q5 + 355 - 2W^r^ [2Q5 + 55 + 2^5]"^ 
215| + - I255WI - 4W| + 2(75^ + m^)Q^S^] [2Q^ + 155^ + 6^5]"^ 

X [2Q5 + 35^ + 3W^r^ [2Q5 + 355 - 2^5]"^ [2Q5 + 55 + 2^5]"^ 
^ [(355 - 2W5)(55 + 2W^){S^ + W5) + 2(5^ + m^)Q^S^] 



6 



2Q5 



[2Q5 + 355 + 3^5]^^ [2Q5 + 355 - 2^5]"^ [2Q5 +s^+ 2w^r^ . 



(6.11) 



We have checked that our full arbitrary group solution correctly reproduces the one loop propaga- 
tor corrections of (13. 6p . Indeed the gluon propagator remains suppressed in the zero momentum 
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limit at one loop. Aside from indicating how involved the final expression for the propagator 



is, the main point of (j6.1ip is to illustrate which of the form factors enhance. From (j6.ip the 
initial terms of the small momentum expansion of represents the gap equation. Therefore, in 
this limit when the gap equation is realised is effectively O ((p^)^) which corresponds to the 



same situation with the Faddeev-Popov ghost, a;° and pj^ fields giving rise to the enhancement 
of these fields. From (j6.1ip the situation is similar since several of the amplitudes have an overall 
factor of Q^. Hence these colour channels will enhance whilst the others will not in the zero 
momentum limit. More specifically the leading order behaviour of the transverse part of both 
bosonic ghost propagators in the — )• limit for an arbitrary group is 



{e^{p)c\-p)) 



{pfip)pfi-p)) 



4^2 



+ 



872 



jabe jcde 



pAp) 



(6.12) 



So one colour channel in addition to those of the original propagator enhances. Although 
we have used the numerical values deriving from the one loop corrections to this enhanced 
behaviour is more general. If instead we examine the leading behaviour of ()6.1ip but for a 
general colour group then we find 
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(6.13) 



with the other form factors being non-singular in the Laurent expansion in Q^. This implies 
that the leading behaviour of the transverse part of the ^'^^ propagator is 



(C(p)C(-p)> 



^ac^bd ^ad^bc j;abe rede 
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(6.14) 



So, for instance, when the gap equation is realised in the three dimensional Gribov-Zwanziger 
Lagrangian, then the corresponding propagator will also be enhanced together with the 
Grassmann fields with the same colour group structure as four dimensions, (I6.12p . Interestingly 
for ^^'^ the colour channel which dominates in the infrared is that which is antisymmetric in 
the colour indices of the field itself. In other words the colour projected field /"^'^C^'^ does not 
enhance at one loop given the relative coefficients of ()6.14p . 






(T 



Figure 4: Several two loop topologies with exchange. 

Having established that and p'^ both enhance we can consider the implications of this 
for the static potential. These are better candidates for a confinement mechanism since they 
are bosonic fields and therefore can be exchanged between coloured static sources. However, the 
same points that were noted in [50] equally apply here in that and /?°^ do not couple directly 
to a (static) field such as a quark or gluon. In the context of the static potential, the colour 
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source only couples to the gluon field, A^, in (|2.3|) . However, one can have a single exchange 
indirectly, for example, through various topologies such as those illustrated in Figure 4. The 
first graph involves source ^"'^ vertex corrections whereas the second graph is a correction to 
the mixed 2-point function. If we consider the first in an extended computation one can replace 
the single by the enhanced propagator whilst the second graph could be ignored since that 
is part of our earlier matrix inversion. Given (j6.12p this exchange will actually give an 0{a?) 
contribution to the static potential and not an 0{o?) one. Although this is clearly contrary to 
a perturbative approach it is a qualitative indication that a reordering of the series will emerge 
even in a Dyson Schwinger approach as already discussed in [50j. Moreover, a similar reordering 
to produce the renormalization group invariant effective coupling constant freezing occured in 
(j3.14p . In the context of freezing it is worth noting that as the Lagrangian channel of both S^'^ 
and p'^ have enhancement, one could define an effective coupling constant based on either the 
gluon vertex or gluon p^^ vertex rather than the ghost gluon vertex which has a freezing 
value the same as (j5.15p in either the MS or MOM schemes. Briefly instead of using either the 
ghost or w^'' 2-point function in the denominator of the definition, one uses instead the first part 
of ()6.ip . However, this is offered as an observation since the renormalization group invariance 
of such a construction is clearly not established and the gluon vertex is certainly not on a 
par with the ghost gluon vertex in the context of non-renormalizability. 

If we take the point of view that our reasoning is at a qualitative level, we can then examine 
some possible implications of enhancement on the static potential. We have computed the 
eight Feynman diagrams contributing to the topology illustrated in the left hand diagram of 
Figure 4 for the static potential formalism. With the propagator of (I2.2ip we find in the 

— )• limit that 
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where, since we are working at two loops, we have put the sources in the adjoint representation. 
Useful in handling the group theory aspects of this calculation was the color .h Form routine 
from [STj based on [75]. However, with the enhanced propagator of (|6.12p we find in the same 
limit 
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+ O (P 



2x2 



where we record the SU (3) result given that the arbitrary group expression would be too cum- 
bersome. The reason for this is that this leading order term does not derive from the enhanced 
part of (I6.12P but instead the 0{l/p^) correction. In other words the dominant part of the 
exchange diagram in the zero momentum limit does not depend on the enhancement. This is 
because the one loop vertex subgraphs corrections are proportional to a structure function. More 
specifically the colour dependence of the vertex is f°-'"^T'^ where the group generator carries the 
indices of the source legs. Closing to form the Wilson loop leads to a trace over these indices 
when the other vertex is included. Since the colour channel of the enhanced part of the 
propagator vanishes identically when the indices of this structure function are contracted then 
that contribution is absent. Indeed we have evaluated the one loop source ^'^ vertex function for 
the static potential momentum configuration exactly and verified that this colour structure is 
the source of the cancellation. Moreover, in this set-up the vertex function in principle involves 
two contributions deriving from the two possible vectors, and p^, which the vertex can be 
decomposed into since it has one free Lorentz index. It transpires that to this order the form 
factor of that of the term is zero leaving only that for v^. Hence, we do not need to consider 
the longitudinal part of the propagator since we are in a static situation where vp = 0. Whilst 
there is the possibility of enhancement from the p'^ propagator through similar topologies to 
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those illustrated on the left in Figure 4, each of the contributing vertex corrections vanish identi- 
cally since the P^i,{k) is contracted with where k is the loop momentum. This occurs because 
the gluon p"^ vertex, which we have not neglected, is essentially proportional to the Landau 
gauge fixing condition. So in fact ()6.16p represents the inclusion of both enhanced propagators 
of (I6.12p . Therefore, we appear to be forced to conclude that the single ^'^ exchange process 
of Figure 4 cannot lead to the dipole behaviour underpinning the linearly rising potential if the 
one loop enhancement of (I6.12P is accepted. 

Indeed considering the higher loop corrections to the source vertex would not appear to 
remedy the situation. For instance, one would have to have a colour structure for the vertex 
which does not involve f°-^'^T^ if the colour structure of (16.12P was preserved beyond one loop. 
Another candidate would be d'^^'^T'^ where d°'^'^ is the totally symmetric rank three tensor but 
it clearly gives zero when contracted with ()6.12p . Whilst one might contrive something to 
circumvent this then the coefficient of whatever this colour tensor is would have to be constant 
in the zero momentum limit. It is not clear at which loop order this could emerge. Alternatively 
when the higher loop corrections to the propagators are computed then it may be the case that 
other colour channels aside from are enhanced. In such a case the colour tensor of (j6.12p 
should be different allowing the enhanced piece to dominate. Clearly such considerations are 
beyond the scope of the current article but suggest that using resummation methods such as 
the Schwinger Dyson technique could probe this in more detail. However, it is worth noting the 
situation with regard to the original observation of bosonic ghost enhancement of Zwanziger, 
|50j . In |50j an over-enhancement was obtained for various colour channels. In order to obtain 
a linearly rising potential it was argued, |50], that the quark gluon vertex develops a Lorentz 
tensor coupling involving a^'^ = Yi^,^'^] together with one momentum vector to reduce the 
overall exchange to a dipole. Since we have considered in essence static gluons we do not have 
the same freedom in Lorentz space to accommodate a tensor coupling. However, provided 
the colour structure did not make the problem trivial then an over-enhancement would clearly 
reduce the power of momentum in (j6.16p at least by another power. Indeed given this it could 
be the case that the higher loop corrections to ()6.ip might also lead to an over-enhanced 
propagator deriving from colour channels other than the propagator one with another colour 
structure to (16.12p . Though we do note that p'^^ enhancement is preserved at two loops in MS 
using (j5.6p similar to c"' and oj'^. To verify this enhancement we evaluated the 212 two loop 
Feynman diagrams of the p^^ 2-point function in the zero momentum limit using the vacuum 
bubble expansion. We recall one advantage of the enhancement feature is that it avoids the 
proliferation of higher order powers in l/p^ which could emerge if one calculated the static 
potential order by order in perturbation theory as indicated earlier. The use of the gap equation 
in being central to this appears to be unavoidably essential to any analysis of studying the 
zero momentum limit. If anything these remarks might only serve to indicate how delicate 
it is to determine the zero momentum behaviour of the localizing bosonic ghost propagators 
via the inversion of the matrix of 2-point functions. Having said this it is important to state 
that we are not ruling out the existence of a linearly rising potential in the Gribov- Zwanziger 
Lagrangian. There are other avenues one could consider aside from the single exchange of Figure 
4. For instance, the enhanced .^^^ propagator can appear inside loop diagrams and hence one 
would have to use a Schwinger Dyson style of analysis to study the implications for the overall 
topologies in the infrared. The enhanced propagator would only be significant at low virtual 
loop momenta. Alternatively the dominant topologies in the infrared could be something such 
as ladder graphs instead of the simple single graphs analysed here. Again that would require 
techniques beyond those discussed here. 
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7 Power corrections. 



Next, we return to our earlier comments concerning the appearance of power type corrections 
of the form ^'^ /p^ at higher loops in the static potential. However, we restrict them to our one 
loop static potential. First, one widely used tool to probe towards the infrared in QCD is the 
operator product expansion. In essence it provides a tool to include corrections to perturbative 
expressions where the corrections involve the vacuum expectation values of gauge invariant 
operators, such as G^^G"'^'^ . In the conventional perturbative vacuum the expectation value of 
such operators is zero but in the true non-perturbative vacuum they acquire a non-zero value. 
Therefore such dimensionful quantities provide the mass scale required to have an expansion in 
inverse powers of the key momentum in the operator product expansion of the particular Green's 
function of interest. One situation where this formalism is applied is to the static potential. See, 
for example, [791 180j. For instance, in [79] the zero momentum exchange of a single gluon with 
a one loop self-energy correction in gluon scattering is examined and a l/{p^)'^ correction is 
produced where the dimensionality is made consistent by the presence of the dimension four 
quantity (G^^.G"'^'^). However, clearly there is nothing to prevent higher order powers of l/p^ 
appearing in such analyses and so it is worthwhile considering (j4.7p in a similar power series 
expansion. Therefore, from ()4.7p we have 
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(7.1) 



Clearly the leading term in this ^"^ /'p^ expansion is the perturbative result of [36^ [371 I38j . 
Moreover, at leading order in the perturbative expansion the next to leading power correction 
is 0(7^/(p^)^) which follows trivially since the full term is a function of 7^ and not 7^ which 
the one loop correction clearly is. Therefore, the leading tree part of the potential mimicks the 
correction observed in |79] although clearly here the quantity used to ensure the dimensions 
balance is 7^ and not Though it ought to be stressed here as was emphasised 

in [79], that this is a short distance approximation to the potential which has not fully been 
accounted for on the lattice. The reasoning is that the confinement force is in principle accessible 
at both low and higher energy scales. A recent exposition on this point has been provided in 
[81j . Though we believe one needs to be careful in this power correction approximation since 
mathematically a short distance expansion of the full potential can only give an insight into the 
r dependence for a limited range of r. So before considering the loop correction, the tree part 
of the potential considered as a power series in 7 gives an interesting insight into applying the 
Fourier transform to coordinate space. As the exact transform is known in ()4.15p . it can be 
expanded in powers of which is also the combination of variables which appear and not their 
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square or higher powers. We find 
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The absence of a term hnear in r is consistent with there being no dipole in the momentum 
space tree potential. However, the one-to-one power series matching clearly breaks down with 
the appearance of a quadratic correction which would ordinarily be associated with a momentum 
space term of l/(p'^)^/^ on dimensional grounds. 



In the one loop correction of (17. ip one observes that a dipole correction appears at next to 
leading order in the power expansion rather than the triple pole at leading order. Parenthetically 
recalling our comments concerning the sign of 7^, then this term would then lead to an effective 
string tension, cr^^, of 
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in the conventional definition of the potential. However, as is evident from our full expression 
there is actually no net explicit dipole term in (j4.7p whose small momentum or large distance 
potential effectively becomes Coulomb-like. If there was a pure dipole present in addition to 
the remaining 7 dependent part, then performing a power series expansion of the full expression 
would mean it would be difficult to isolate its contribution uniquely. Moreover, it would be 
difficult to interpret such a correction in an operator power series context since it would involve 
the square root of the gluon condensate or the vacuum expectation value of a dimension two 
object. In the Gribov-Zwanziger context such an operator could be f"-bcj^afi^ab from the 

equation of motion 
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would effectively equate to the presence of the Gribov horizon condition operator. Though one 
other application of ()4.7p could be to use it as a testbed for examining renormalon style correc- 
tions, like [8OJ, since the gap equation solution (|5.3p is clearly non-perturbative. For instance, in 
the F-scheme we would have a power correction beyond the perturbative contribution of (|4.10p 
in (I4.19P which is 
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where, [36l Ell [38] , 
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A formally similar power correction was observed in the effective coupling of (j5.13p in |,33J where 
the correction had the same sign. In order to compare and for completeness we record that the 
power correction to this renormalization group invariant coupling in the same notation as (j7.6p . 
[331. is 



aefl(p) 



a 



pert 
efr 



(P) 



16p2 




(7.7) 



35 



where 



1 + 



169 11 

"36" ~ y 



Ca + 



3 



PI 



20 



a(^) 



(7.8) 







Figure 5: Leading contribution to operator correlation function. 

Next, we comment on one aspect concerning potentials and that is whether it is possible to 
observe stable bound states. If there were such states then these could possibly be identified 
with glueballs. However, since we have demonstrated that the static potential does not have 
a linear confining part at one loop but a form not dissimilar to that of Vo(r), we expect any 
states formed in, say, a Schrodinger equation analysis of (|4.7p to be unstable. For instance, in 
|67j a confining piece had to be included in that glueball analysis which assumed a gluon mass. 
For a similar calculation in the Gribov-Zwanziger context one would have to implement the gap 
equation for the Gribov mass in some fashion in addition. Aside from this one way of possibly 
quantifying properties of bound state masses is by considering the correlation of operators with 
the same quantum numbers as the bound states along the lines of [13] . For glueballs the obvious 
candidate operators are those involving the field strength given by 

Os = - -^G^C^^- , Oi^" = - \ {^-!^G%G-'^P - G-^'^Gl,^ (7.9) 

since they are gluonic and gauge invariant, where the subscripts 5 and T denote Lorentz scalar 
and tensor operators respectively and the second operator is symmetric and traceless. It is based 
on a similar operator considered in [82j. Indeed Zwanziger has examined the first operator in 
the context of (j2.19p by considering the leading term of the correlator of Os, [Ej. In [13] it was 
suggested that there was evidence for a state with mass squared 2\/Ca7^, o^i" conventions, 
by rewriting the correlation function in a spectral representation using tools such as Schwinger 
and Feynman parameters and searching for physical cuts. However, a full explicit expression for 
the correlator as a function of the momentum and 7 at leading order was not given. Therefore, 
to partly address this we have computed the leading order term of the correlator of 0$ exactly. 
The Feynman diagram is illustrated in Figure 5 where the crossed circled denotes the location 
of the operator insertion with momentum p flowing through it. Defining 



ns(/) = (4vr)2i / cfxe'P^iS)\Os{x)Os{m) 



(7.10) 



we have 



^syP ) = — — - + — tan 



4e 

TrCy 7 , P 



4p2 



p2 



4p2 



+ ^V4C^7^-(p')'tan-i 



^aGai^ - {p'Y 



36 



+ 



2 

1 

4 



In 



3 3^/2 2, 



Ca7^ + 



37r ^/2 2^ 



^2 



V2 
32 



A^A + 0{a) 



(7.11) 



where the divergent contact terms are included for completeness. These are absorbed by the 
canonical contact renormalization for operator correlation functions but we note that as there is 
a Gribov parameter present this contact renormalization actually has a mixing aspect which is 
not unexpected. Given this one can regard the finite part as the leading piece of the correlation 
function since the first appearance of an explicit coupling constant is at next order, ignoring the 
coupling constant implicit in 7, (15. 3p . We have checked that there are no poles in the expression 




and hence regard this correlation function 



or 4a 



at obvious places such as yfCjci^ , 
to be regular as a function of . However, one can see that there is a cut at the same value as 
observed in [l3j which is ]? = 2\JCa^'^ and this is at the same point as the cut in ()4.7p . Equally 
it is elementary to verify Zwanziger's other observation in |13j that there are unphysical cuts at 
p2 = lb 4i-v/CA7^ in our conventions similar to (14.70 . One final point concerning the structure 
of n5'(p2) in relation to both glueball states and the physical cut structure and that is that 
the same conclusion would be obtained if other operators with similar properties to O5 were 
considered. For instance, the correlation functions of both the operators Tr \(J)p,Gya) {D^^C^")] 
and Tr [{D^DyG^p) {D^'D''G''P)] with themselves, where Gf,y = G^^T", will produce the same 
physical and unphysical cuts as (j7.1ip . These Lorentz scalar operators are both gauge invariant 
and have the same number of leading gluon legs. 
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Figure 6: Higher order graphs contributing to operator correlation functions. 

For higher loop graphs contributing to the gluonic operator correlation functions it is straight- 
forward to determine the location of the physical cuts. Such graphs are illustrated in Figure 6 
where the dots denote more and more gluon propagators. When there are an odd number of 
gluons then there are unphysical cuts but no physical ones. For an even number of gluons, say 
2n, then there are both types of cuts with the physical ones being at = 2v?\/Ca'^'^ ■ The 
pattern of unphysical cuts cannot be written in as compact a formula. As noted by Zwanziger, 
|13j . for the two gluon case the cuts additional to the physical ones are at = it M^fGATi^ . It is 
straightforward to record those for the lowest order cases. For the three gluon situation the only 
cuts are at = zb 1-^/6^7^, = ± ^iyfCAri^ and = (it 4 it 3i)\/C^7^ where all possible sign 
combinations are taken in the final off axis cut. For four gluons, the unphysical cuts are now at 
= lb M^fCAri^ , = ± Voi^fCAri^ and = 2{± 4 it 3i)^/CA'y'^■ For other Green's functions, 
such as the higher loop corrections to (|4.7p . we expect the physical (and unphysical) cut struc- 
ture to be the same. Returning to Figure 6 and, for instance, ignoring the presence of internal 
vertices for the moment, the next physical cut in (17.1ip will be at = 8-v/Ca7^. If instead of 
(j7.1ip one considers the correlation function of a gauge invariant operator involving three gluons 
in the leading leg term, such as /"'^'^G'^j^^G^'^^G'^^'^ which was considered in [83], then the first 
cut would be at = 8\/GA'y'^ with the subsequent one at = 18y/ Ca'j'^- In [83] QCD sum 
rules were used to estimate the masses of various glueball states. For instance, the two gluon 
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scalar glueball mass was M2g = (1.50 it 0.19)GeV and that for the three gluon bound state was 
M^g = S.lGeV with no errors quoted, [83]. Intriguingly the mass ratio is M2,g/M2g ~ 2. If we 
take the ratio of the leading physical cuts for the appropriate gauge invariant gluonic operator 
correlation functions then we have the same ratio deriving from the expression of the ratio in the 
crude form ^^8/2. By contrast, however, if one was working with gluons which had an ordinary 
massive propagator instead of the Gribov width, then one would have a canonical cut structure 
in the corresponding Feynman diagrams. Overlooking the obvious loss of gauge invariance for 
the moment, then the corresponding ratio for M^g/M2g would be 3/2. This seems to be too 
far away to accommodate the mass ratio of [83]. Returning to the Gribov-Zwanziger case, an 
estimate was also given in [83] for the tensor 2"'""'" state which was M2++ = (2.0 ± 0.1)GeV. 
This gives the ratio M2++ /M2g ~ 4/3. Such a ratio can be accommodated in terms of ratios of 
physical cuts, such as \/32/18 in crude form. Similarly, M^g / M2++ can be accommodated by 
the ratio -^18/8. However, we acknowledge that the justification of these latter ratios from a 
simple diagram argument similar to the earlier one seems to be tenuous. This is partly because 
(j4.7p can only really be regarded as being relevant for spinless zero angular momentum bound 
states such as the lowest lying two and three gluon states. For states such as 2"'""'" one requires 
more detailed structure beyond the simple radial dependence provided in (|4.7p . such as that 
considered in depth in [67]. It would be interesting, though, to see if the ratio of 4/3 could 
be justified from that point of view in an extended potential incorporating spin and angular 
momentum. However, even though we are primarily concentrating on Yang-Mills here the issue 
of degeneracy will arise at some point for the higher states. Whilst not as involved a problem as 
when quarks are present, which would clearly extend the number of potential bound states, re- 
solving any mixing will require a much more detailed analysis of a potential derived from ()2.19p . 
Another test of the ratio hypothesis would be the mass of a four gluon state which ought to be 
roughly three times that of the lowest state if our reasoning is sound. However, such a state 
has not been extensively studied, for example on the lattice, as far as we are aware. Finally, in 
this simple cut analysis we need to temper our remarks with the fact that with the inclusion of 
higher loops it is not inconceivable that the cut locations will be shifted by corrections. So our 
method of estimating mass ratios should only regarded as a rough guide. 

Whilst trying to understand the relative sizes of various glueball states is an interesting 
exercise in itself, ultimately one has to eventually estimate one of the masses which is a non- 
trivial task. We draw attention to several possibilities in the context of (12.19p . In each case one 
has to somehow fix a mass scale in relation to a measured quantity from other methods such 
as the lattice in order to obtain a numerical value for yfCjci^ which is the fundamental mass 
parameter in (|2.19p . One approach to extract a glueball mass would be to consider the power 
corrections to (|7.1ip and apply sum rule technology akin to that used in |82l [83] where the gluon 
condensate is the quantity used to fix a scale. In [82] a tachyonic gluon mass was introduced to 
account for apparent discrepancies with experimental data. Performing the detailed analyses for 
both these approaches is clearly outside the scope of the current article. For example, a proper 
sum rule analysis would first require a reworking of the original operator product expansion 
but using the Lagrangian of (|2.19p rather than the usual QCD Lagrangian. However, we have 
computed the leading power correction to (jT.lip in the context of the gluonia channels of VLsiff'^ 
considered in |82j. Therefore, using the same moment definition and notation as |82j we have 



n5(M^) = (parton model) 1 



M4 



) 



(7.12) 



where 




(7.13) 
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and = /n is finite in the limit of large and n. As a check on our derivation of this power 
correction using FORM, [57], we have replaced the gluon propagator of (|2.19p by the tachyonic 
gluon mass propagator of [82] and correctly reproduced the power correction quadratic in this 
mass recorded in [82]. As we are dealing with a gauge invariant operator correlation function in 
the Gribov-Zwanziger case the first power correction is quartic with the correction numerator 
being related on dimensional grounds to the gluon condensate. Although in the static potential 
case the leading order function was a function of 7^, the loop correction produced 0(7^) terms in 
the power series expansion. Such a scenario could emerge for Jisip^) when the 0{a) correction is 
computed. Though this is currently beyond the scope of this article as well as using sum rules to 
estimate scalar or tensor glueball masses akin to |83] . Finally, for the correlation function of the 
tensor operator Ojf, similar power corrections should emerge in each of the scalar amplitudes of 
its Lorentz decomposition. We record the full expressions for the tensor case for completeness, 
and for comparison with [82], in appendix C as they are similar to (|7.1ip . However, there is 
nothing formally different from the scalar case above having, for example, the same physical 
cut. 

Another approach would be to make direct contact with a non-zero vacuum expectation 
value, whose numerical value is known, by explicit computation using ()2.19p itself. This is 
possible since now the gluon propagator is not massless, (|2.2ip . and so vacuum expectation 
values of gluonic operators are non-zero. Such an approach has already been used in [SH [85l |86] 
but where the gluon is assumed to have a canonical mass term. Ignoring the first of these three 
papers, since it appears to have a result inconsistent with the latter two, an estimate for the 
gluon mass was deduced by comparing the leading order gluon condensate value with the quark 
current correlator in the high energy limit at the same order. However, we can also consider 
the dimension two condensate based on the operator ^(Ap^ which has been measured on the 
lattice in [721 ES] • At leading order we have 



So, for instance, at leading order we have the formal relation for a glueball mass, if we regard 
our lowest cut as a glueball mass, of 2VCU7^ = (167r/3)(i(^;i)2). Taking the lattice estimate 
of {{A^f) = 3.1GeV^ from [7^, for example, then this rough way of estimating would give an 
unrealistic glueball mass of S.lGeV. Aside from this one ought not to overlook the gap equation, 
(j5.2p . which is assumed to be valid at this order of approximation. If one uses representative 
numerical values for \/Ca7^ derived in this way and a reasonable estimate for Aj^ of, say, 
^MS ~ 300MeV then one would obtain a negative value of which is clearly unacceptable. 
Though this is partly due to the Landau pole problem of the running coupling constant. (We 
find similar conclusions when the gluon condensate is used similar to [85].) Of course, this simple 
exercise has been recorded to merely illustrate some of the potential difficulties in estimating a 
value for the underlying mass parameter of p.lOp . Indeed we have only considered leading order 
and including higher order corrections will in principle change the estimates. However, the full 
vacuum expectation value should have a non-perturbative piece which would need to be properly 
incorporated into any deeper analysis as well as dealing with infrared issues which are known to 
exist in the case of the gluon condensate. The lack of consistency with the Gribov gap equation 
to the order we considered, clearly indicates that the approximations assumed here could not 
be validated in a consistent way. By contrast there is no equivalent gap equation constraint on 
the fundamental gluon mass used in the analyses of [HH [851 ES]- Also, it may be the case that 
the Gribov gap equation would require a more complete function of a beyond the perturbative 
approximation. Therefore, whilst carrying out this rough leading order analysis has led to a null 
conclusion, we feel it is useful to include it as a moderating point of view since it does perhaps 
illustrate the difficulty in producing a reasonable estimate for the underlying mass parameter of 
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()2.19p . Indeed it may also be indicative that the use of condensates to fix a mass scale may not 
be compatible when a Gribov mass is present. 

8 Discussion. 

The main result of this article is the explicit construction of the one loop static potential using 
the Gribov-Zwanziger Lagrangian. The original motivation was to study a gauge invariant object 
related to the non-perturbative structure of a formulation of a non-abelian gauge theory which 
has properties suggesting it describes a confined gluon. Indeed one aim was to see what functions 
of momentum appeared in the explicit final expression. Whilst a linear potential clearly does 
not emerge in the final result, there is an intriguing hint with the presence of a pure dipole term 
but which has a compensating term in the zero momentum limit for positive 7^. Ultimately we 
have to conclude that the usefulness of this approach to the static potential is that it will in 
principle closely match the full potential if one adds successive loop corrections to (14. 7p . So it 
would appear that in the present context a linearly rising potential for a significant range of r 
arises out of a truly non-perturbative effect. Indeed it is hard to see how the compensating term 
could be split from the pure dipole in higher loop corrections to leave a net dipole as well as 
no net higher order momentum poles. Therefore, some non-perturbative feature must be taken 
into account. 

In the context of this field theoretic approach in the Gribov set-up the most promising 
candidate for this is the (non-perturbative) Gribov gap equation. To this end we revisited the 
behaviour of the bosonic p'^ and localizing ghosts which dominate in the infrared and are 
responsible in essence for implementing the horizon condition in (I2.19p . We have been able to 
derive enhancement for both bosonic localizing ghosts, together with the colour structure but 
it differs from that derived in Zwanziger's Schwinger Dyson analysis, [50]. Moreover, we do not 
find the over-enhancement discovered in [50j . We have discussed various possibilities that might 
lead to similar behaviour but this would at least require a higher loop computation. One lesson 
appears to be that the inversion of the matrix of 2-point functions to obtain the propagators 
is intricately tied to the colour group structure of the corrections. For instance, to extend the 
perturbative analysis to the two loop level is in principle possible in the zero momentum limit 
but requires the computation of over 1000 Feynman diagrams since one has to consider the full 
matrix of 2-point functions in the sector. This would at least give some indication of 

the zero momentum behaviour of the colour channels not present in the original Lagrangian and 
whether the gap equation emerges in the corrections similar to (16. ip . It could be the case that 
additional enhancement arises in other channels or an over-enhancement akin to that of |50] . 
The aim of such an analysis would be to see if a dipole dominated the exchange between the 
static colour sources of the formalism but in such a way that higher order poles were excluded. 

Further, one underlying feature of the static potential formalism is the absence of a direct 
coupling of the source to either p'^ or which prevents the direct sing le pf or field exchange 
with an enhanced propagator being the simple explanation for a linear potential, provided such 
a vertex did not depend on the exchange momentum. Though one would naively believe has 
to play some role in this way since from its equation of motion it is a non-local projection of the 
gluon. Perhaps the original Wilson loop static potential formalism would need to be reconsidered 
in the Gribov-Zwanziger case due to the restriction of the path integral to the Gribov region 
of configuration space. From another point of view, although the Gribov-Zwanziger Lagrangian 
is consistent with the Kugo-Ojima confinement criterion, [341 135j . that is a necessary but not 
sufficient condition for confinement which would imply that an additional feature might be 
necessary in (|2.19p to obtain a confining potential. For instance, the Gribov construction is 
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founded on the infinitesimal behaviour of the gauge fixing condition, [10], and the less local 
aspects of that construction might need to be included now. 

Throughout the article we have concentrated on what is referred to now as the conformal or 
scaling solution. In recent years there has been interest in an alternative point of view which is 
called the decoupling solution, [22l[23l[2ll[25l[26l[271|28l[29]. Essentially this scenario differs from 
the gluon suppression and ghost enhancement properties of the conformal solution, in having no 
Faddeev-Popov ghost enhancement and the gluon propagator freezes to a finite non-zero value 
excluding suppression. The evidence for this behaviour derives from both lattice gauge theory 
and Schwinger Dyson analyses, [22l [231 [Ml [Ml EZl [28l [29] . As yet there is no definitive 
concensus as to which of the conformal or decoupling solutions is the correct picture of infrared 
Landau gauge Yang-Mills theory. 

One proposal to explain the lack of suppression and enhancement in the Gribov-Zwanziger 
Lagrangian is the condensation of the BRST invariant operator 

^ab^abf, _ ^;^ab^abf^^ [SSI [53], 

using the notation of the original localizing fields. With a non-zero value for the vacuum ex- 
pectation value of this operator the observed behaviour of the gluon and Faddeev-Popov ghost 
propagators could be accommodated. However, to include such an operator one applies the 
local composite operator formalism developed in [871 [881 [89] . Briefly the method introduces a 
new source coupled to the dimension two operator and then constructs the effective potential 
of the operator. Studying the minima of this effective potential one observes that the perturba- 
tive vacuum solution is unstable in favour of a vacuum solution where the operator condenses. 
Hence a new dynamically generated mass is introduced which modifies the propagators of the 
appropriate fields. In the application to the Gribov-Zwanziger Lagrangian, [52l[53], both sets of 
localizing ghosts acquire a mass and induce extra mass dependence in the gluon and Faddeev- 
Popov ghost propagators. Moreover, the corresponding gap equation is insufficient to produce 
a Faddeev-Popov enhancement. 

This raises several points for our static potential calculation. First, we have taken the point 
of view that the p^^, ^"'^ and co^ localizing fields are purely internal fields with no coupling to 
external sources in the way the gluon does. Therefore, it seems unclear how to incorporate all 
the dynamically modified propagators at the outset since including them would require coupling 
an operator of internal fields to an external source. Moreover, to have a homogeneous renormal- 
ization group interpretation in the local composite operator formalism one has to allow for the 
generation of the square of this external source. However, even if one relaxed such assumptions 
or used the modified propagators directly in the analogous computation of (|4.7p then it seems 
difficult to see how a linearly rising term could emerge in the corresponding static potential due 
to the absence of enhancement. 

This leads on to the other point of view we examined and that was the enhancement of 
p'J^ and It now appears evident, at least in the conformal solution, that the enhancement 
of the Faddeev-Popov ghost, p'^, .^^^ and u)'^ are on the same footing and inextricably linked. 
Therefore, if there is a loss of enhancement in the decoupling solution for the Faddeev-Popov 
ghost it would seem inevitable that this would be the case for the localizing fields too. Whilst 
we analysed the simple two loop contribution of Figure 4 to the static potential with enhanced 

and pI^^ propagators, it was insufficient with the enhancement to produce a dipole behaviour 
at one loop primarily due to group theory considerations. If there was an non-enhanced 
propagator in the decoupling case then there would appear to be no emergence of anything 
like a dipole behaviour in the zero momentum limit. However, the behaviour of the and 
^'^ propagators in the infrared for the decoupling solution has not been established yet. If the 
decoupling solution is established as the correct description of the infrared properties of Landau 
gauge Yang-Mills it would be interesting to see how the linearly rising potential emerges in 
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the field theory context. For instance, in [19] the Wilson loop was considered in the context 
of (j2.19p and it was shown that the string tension of a linearly rising potential depended on 
the zero momentum value of a particular combination of the gluon and Faddeev-Popov ghost 
propagator form factors. Briefly, if the gluon propagator froze to a zero or non-zero value but 
the Faddeev-Popov ghost was enhanced then a linear potential emerged. If the argument of [19j 
remains valid in the decoupling scenario then the lack of Faddeev-Popov enhancement would 
seem to exclude a linear potential in that case. 
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A Group theory. 

In this appendix we discuss the decomposition of products of group generators into the basis of 
Casimirs. For a comprehensive review of such Casimirs we refer the reader to [78j. In section 6 
we introduced the two sets of coefficients {aj} and {bi} where 
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and 
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ade rbce 



jabcd 



japejbqe^cdpq ^ ^^^ab^cd j_ (^^ac ^bd ^ ^ad^fec^ ^ (^jace jbde ^ jade jbce\^ ^ ^^^abcd _ (^ 2) 

These coefficients can be computed by the projection method. By this we mean that each tensor 
on the right hand side is used in sequence to multiply each equation. Using the properties of 
the Lie algebra this leads to a set of linear equations for the Oj and hi which can be solved by 
simple matrix inversion. It transpires that the matrix which needs to be inverted is the same 
for both ([XT]) and ([A2]) and is 



M 



2Na 

2Na 2Na{Na + 1) 
2CaNa - 2CaNa 

V IcIna IcIna 



2CaNa 

c 

^CINa 



2CaNa IC\Na 


^abcd^abcd 







-eClNA \ 
'CINa 




To determine and bi the inverse of M multiplies the respective vectors 
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Hence, we find 
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a2 



1U{11Na - 8) - 1080CINa{Na - ^)df'''^(f^P''d'f"' 



abed jabcd 



03 



04 



108Na{Na - 3)[12(iVA + 2)d^l^^'d^l^'^ - 25^47X4] 
[12{Na + 2)d"^'"^df'"^ - 25C\Na] 

54CAiy4(iy4 - 3) 

[216(iVA + 2)(i^''"^4^P^4''P^ - 125C%Na - 360CK"'^<"'^] 
18[12(Ar4 + 2)d^J^3^d^J^^ - 2bC\NA] 



(A.5) 



and 
bi = 
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These expressions have been derived by making use of the color . h package of j571 [75] . Both sets 
of coefficients can be evaluated exphcitly for SU{Nc) and for completeness we note the values 
of the various Casimirs in this instance are 
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where the last expression is given in |9U| . The set {aj} also require d'^^'^d^^^'^d"^^'^ which we have 
evaluated directly using the SU (Nc) identity of fOT] 
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and relations for products of the structure functions and the totally symmetric rank 3 tensor 
^afec^ [91]. As an intermediate step we have 
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for SU{Nc) which we have checked correctly reproduces the final expression of dAJl). Hence, we 
found 

dt'df'dt' ^ N^ [N^ + 135iV,^ + 324] 

Na 216 ■ ■ ) 

We have checked the consistency of this expression by explicitly evaluating the left hand side 
for both colour groups SU{2) and S'f/(3). To do this we used the explicit respective 3x3 and 
8x8 matrix representations of the adjoint group generators using Form. Although it is not 
required here, as a corollary we have also determined the SU{Nc) value for d"^'^'^'^^ d"^'^'^^^ which 
is introduced in [Ml- We find 
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Equipped with the SU{Nc) values of these Casimirs we have 
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It is worth noting that some of the numerator factors as well as the denominator factors of Oj, 
which are Na, (iYj - 3) and [12{Na + 2)df'"^df'''^ - 2bC\NAl have zeroes at Nc = l,2 and 3. 
Therefore, one cannot directly evaluate for the latter two values of but must derive the 
Nc dependent expressions first before determining numerical values. 



B Transverse parts. 



In this appendix we collect the explicit one loop MS expressions for the transverse parts of the 
2-point functions of (13. We have 
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Taking the zero momentum limit of each expression we find 
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C Longitudinal parts. 



In this appendix we provide the one loop MS expressions for the longitudinal parts of the 2-point 
functions of (13.11). We find 
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As in the previous appendix, taking the zero momentum hmit of each expression we find 
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Interestingly the 0(7^) part of is equivalent to the 0(7^) term of X in the same limit. 



D Tensor operator correlation function. 



In this appendix we record the explicit form of the correlation function of the Lorentz tensor 
operator Oj!^ which is defined by 



(D.l) 



We decompose this into a similar basis of Lorentz tensors to that used in [8^ with the only 
difference being that we work completely in d-dimensions and not four dimensions, as we use 
dimensional regularization. We have 
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To determine the scalar amplitudes Hf{p'^) we use a projection method. This involves multi- 
plying the original correlator by a Lorentz tensor which is a linear combination of the three 
basis Lorentz tensors defining the decomposition, ()D.2p . The construction of each of the three 
projection tensors is similar to that used in appendix A for the decomposition of colour group 
tensors. We first construct a matrix where the entries are determined by multiplying ()D.2p by 
each basis tensor in turn. This produces the d-dependent matrix 
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The coefficients of the basis tensors in each of the three projectors are then determined from 
the inverse which is 
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Equipped with this we find the three scalar amplitudes are 
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Again we have included the divergent parts of the correlation functions which are absorbed by 
a contact renormalization. We note, though, that there is only mixing for the first amplitude. 
Aside from several terms involving the factors l/(p^)^ and l/(p^)^ the actual functions appearing 
in each of the amplitudes are the same as for 115(1?^). So, for example, the cut structure is the 
same with a physical cut at = 2-^/0^7^ • Repeating the same moment calculation for the 
scalar operator of section 7, we find that 



and there are no O 
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corrections for the other two scalar amplitudes. Further, if one 



considered the correlation function of the energy momentum tensor, as in [82], then the same 
power correction as 11^ (M^) would emerge. 
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